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Abstract. This paper is concerned with the Becker-Doring (BD) system of 
equations and their relationship to the Lifschitz-Slyozov- Wagner (LSW) equa- 
tions. A diffusive version of the LSW equations is derived from the BD equa- 
tions. Existence and uniqueness theorems for this diffusive LSW system are 
proved. The major part of the paper is taken up with proving that solutions 
of the diffusive LSW system converge in the zero diffusion limit to solutions of 
the classical LSW system. In particular, it is shown that the rate of coarsen- 
ing for the diffusive system converges to the rate of coarsening for the classical 
system. 



1. Introduction. 

In this paper we shall be concerned with the Becker-Doring (BD) system of equa- 
tions [5] and their relationship to the Lifschitz-Slyozov- Wagner (LSW) equations 
[TT1 [22] . The BD equations describe the time evolution for a mean field model of a 
population of particle clusters. Particle clusters are characterized by their volume 
I which may take integer values I = 1,2.... The monomers of volume I = 1 play a 
distinctive role since they mediate interactions between particle clusters of different 
volumes. In particular, a cluster with volume £ can become a cluster of volume £ + 1 
by addition of a monomer, or become a cluster of volume £ — 1 by evaporation of a 
monomer. Finally, conservation of mass imposes a constraint which is non-local in 
£. 

The BD equations are determined by the rates at which a cluster of volume £ 
becomes a cluster of volume £ + 1 or £ — 1. If c\(t) is the monomer density at time 
t, then a cluster of volume £ becomes a cluster of volume £ + 1 at rate a£C\(t), 
and a cluster of volume £ evaporates a monomer at rate bi to become a cluster of 
volume I — 1. The net effect of these two processes yields a flux Jg of ^-clusters 
to £ + 1-clustcrs which depends on the density cg(t) of ^-clusters and cg + i(t) of 
£ + 1-clusters, as well as the monomer density Ci(i). In the BD model it is given 
by the formula Ji = agcici — b^+iQ+i. The BD equations for the time evolution of 
the system are therefore given by 

(i.i) ^ = Jt-i-Jt, e>2, 



(i.2) y><(*) 
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Global existence and uniqueness theorems for the BD system (|1.2p were 

proven in the seminal paper of Ball et al pQ under fairly mild assumptions on the 
rates at, be, £ = 1,2..., the main one being that ag should grow at a sub-linear rate 
as £ — > oo. 

It is easy to see by solving Ji = 0, I > 1, that (|1.2p has a family of 

equilibrium solutions parametrized by the monomer density c\ . In order to classify 
the equilibrium solutions one needs detailed knowledge of the asymptotic behavior 
of the rates ae, be as £ — ► oo. In this paper we shall assume that ae, bp are given 
by the formulae 

a e = ai f/ 3 , ai > 0; b e = a e (z s + q/£ 1/3 ), 



(1.3) 



q > 0. 



These values for ae , be may be obtained by deriving the BD equations as mean field 
equations for a 3 dimensional microscopic system [20] . 

With ae, be given by (|1.3p . one sees that for any c\ satisfying < c\ < z s there 
is an equilibrium solution C£ = Qec\, t = 1, 2.... For large £ the coefficient Qe has 
the asymptotic behavior 

1 



(1.4) 



2z, 



oo. 



It follows from (|1.4|) that there is an equilibrium solution corresponding to c\ = z s . 
We shall denote its density fll.2|) by /9 C rit- It was shown in [1] that if p < p cr it, then 
the solution of (jl.ip . (|1.2p converges strongly at large time to the corresponding 
equilibrium solution in the sense that 

oo 

lim Y £\c e (t) ~ Q e c[\ = 0. 



If p > /9 cr it then ce(t) converges weakly to the equilibrium solution with maximum 
density 



(1.5) 



lim ce(t) = Qgz{, £>1. 



Equation (| 1 . 5|) shows that for p > p cr i t particles from the excess density p— p cr it over 
equilibrium concentrate at large time in clusters of increasingly large volume. This 
is the phenomenon of coarsening, and it is an important problem to understand the 
rate at which the volume of clusters from the excess density increase with time. 

A mechanism for determining this was proposed by Penrose in [20j . where he 
argued that the large time behavior of clusters from the excess density is governed 
by the LSW equations, 



(1.6) 



(1.7) 



dc(x, t) 
dt 

xc(x, t)dx 



d_ 

dx 



L(t) 



1/3" 



c(x, t) 



x > 0, 



Pcrit- 



The parameter L(t) in (|1.6p . which turns out to be a measure of the average cluster 
volume, is uniquely determined by the conservation law (|1.7p . A mathematically 
rigorous relationship between solutions of the BD system (jl.ip . (|1.2|) and LSW 
system (|1.6|) . (|1.7p was established in [TD] for the case z s = when p clit = 0, and 
in [T3] for the case z s > when p crlt > 0. Much still remains to be done however, 
in order to understand this in the precise way envisioned by Penrose j2Tj] . 
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In this paper we take a rather different approach from |10[ 113] in attempting to 
establish a relationship between solutions of the BD and LSW systems. Following 
[2"T] . we first observe that the flux Ji can be written as 



Observe from (|1.5j) that the function L(t) in (|1.8|) satisfies L(t) — > oo as t — > oo. 
Denoting by D the forward difference operator and D* its adjoint, then (|1 . ip is 
given by 



where c(£,t) = ci(t) in Evidently if we drop the first term on the RHS of 

(|1.9p we obtain a discrete version of (11.61) . 

We still need to account for the discrepancy in the conservation laws (|1.2|) for 
the BD system and (|1.7|) for the LSW system. Now from (|1.5[) the bulk of the 
equilibrium density p cr ;t is concentrated in clusters of volume 0(1) as t — > oo. 
Further, by choosing £* large enough we may make c{£* ,t) arbitrarily small for 
all large t, and the density concentrated in clusters with volume larger than £* 
arbitrarily close to p — p cr it- It may seem reasonable therefore, that in order to 
understand the large time behavior of the BD system, one is justified in setting 
c(£*,t) = 0, t > T, for some large but fixed I*. This assumption is evidently 
unphysical since it cuts off interaction between clusters of volume less than £* and 
clusters of volume greater than £*. One consequence of it is that the resulting 
conservation of mass for clusters of size greater than £* implies C\{t) > z s for t > T. 
There is so far no rigorous proof that solutions of the BD equations with p > p cr i t 
satisfy ci(t) > z s , t > T, for sufficiently large T. A proof of it would seem to be 
necessary in order to rigorously justify the relationship between the BD and LSW 
systems. In this connection one should also note that in |13j it is necessary to make 
an assumption -equation (1.33) of the paper-in order to prove that a scaling limit 
of the BD system is a solution of the LSW system. Equation (1.33) is however 
not a boundary condition, but essentially an upper bound assumption on the rate 
of coarsening for the BD model. Some interesting heuristic explanations [5] have 
recently been given to describe the transition, from the relaxation to equilibrium 
phase of the BD time evolution, to the coarsening phase governed by the LSW 
equations. These may be helpful towards constructing a fully rigorous theory. 

We shall assume now that there is a time T such that for t >T, the interaction 
between clusters of large volume and clusters of 0(1) volume is negligible, and 
even make the significantly stronger assumption that there exists £* > 1, such that 
c(£*,t) — 0, t>T. We may therefore study the large time evolution of the cluster 
density for particles from the excess density by solving the initial value problem 
for (|1.9p . subject to conservation of mass for particles with volume larger than £*. 
Without loss of generality we may normalize I* = 1, T = 0, whence the problem 
becomes solving (|1.9p subject to the Dirichlet boundary condition c(l, t) = 0, t > 0, 
given initial condition q(0) = 7^, £ > 1, and conservation law 



J e = [b e c e - b e+1 c e+1 } - a iq [l - {£/L(t)} 1/3 ]c e , 




c 1 (t)=z s + q/L(t) 1 / 3 . 



(1.9) 



dc(£, t) 
dt 



D*D[b lC {£, t)] - a ig D* [(1 - {£ / L(t)}^ 3 )c(£, t)], £> 2, 



OO 



(1.10) 
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In §2 we prove global existence and uniqueness of a solution to this problem (The- 
orem 2.2), and obtain some properties of it. 

Next we consider a continuous version of (jl.9|) . Setting x ~ I and replacing finite 
differences in (| 1 .9[) by derivatives, we obtain the equation 



(1.11) 

dc(x,t) d 
dt dx 



e(l + x/e) 1/3 c(x,t) 



d_ 

dx 



x > 0, 



with the parameter e = 1 in p. lip . In p. lip we have replaced the coefficient 
be = ai(q + Zgl 1 ^ 3 ) of (| 1 . 9[) by the diffusion coefficient (1 + x) 1 ' 3 . This is simply for 
convenience, since the only properties of the diffusion coefficient in p. lip which we 
will use are that it behaves like x 1 / 3 for large x, and is bounded away from for 
small x. An equation similar to (jl.lip was obtained in |21j . In §3 we prove global 
existence and uniqueness of solutions to (jl.lip subject to the constraint (|1.7p and 
the Dirichlet boundary condition c(0, i) = 0, t > 0, provided the nonncgativc initial 
data c(x,0), x > 0, satisfies the intcgrability property, 



(1.12) / (1 + x) c(x, 0)dx < oo. 

Jo 

This is the content of Theorem 3.1. In [21| an analogous local existence and unique- 
ness theorem is proved (Theorem 6 of [21]), for which one must also assume the 
point-wise estimate sup x>0 c(x, 0)/(l + x) 2 / 3 < oo on the initial data. 

There are several existence and uniqueness theorems in the literature for the 
LSW system (|1.6p . (|1.7|) and its generalizations. Collet and Gordon [3j have proved 
existence and uniqueness for a system with globally Lipschitz coefficients, under the 
condition that the initial data c(-,0) is integrable and has finite first moment. The 
corresponding result for the LSW system-in which the coefficients are not globally 
Lipschitz- was proved by Laurengot [5] ■ Niethammcr and Pego |15| proved existence 
and uniqueness for the LSW system under the condition that the initial data c(-, 0) 
is a Borel measure with compact support. In |16| they were able to extend their 
method to the case when the initial data c(-,0) is only assumed to be a Borel 
measure with finite first moment. 

In §4, §5 we shall show that the solution of the diffusive LSW system (|1.7p . 
(|l.lip with given initial data, converges as e — > on any finite time interval to the 
solution of the classical LSW system with the same initial data. One hopes that 
these results will shed some light on the problem of understanding the mechanism 
of coarsening for the BD and other models of Ostwald ripening. It has been argued 
in the physics literature Q~2] that in models of Ostwald ripening a diffusive LSW 
equation occurs. Although the effect of the diffusion term vanishes at large time it 
acts as a selection principle, so that the large time behavior of the model is governed 
by the unique infinitely diffcrentiablc self-similar solution of the LSW model. 

The fact that the effect of the diffusion term is expected to vanish at large time 
is closely related to the dilation invariance property of the LSW system. Thus if 
c(x,t), L(t), t > 0, is a solution to (II. 6p . (|1.7p . then for any A > one has that 
c\(x,t), L x (t), t > 0, defined by c x (x,t) = Xc(Xx,Xt), L x {t) = L(Xt)/X, t > 0, 
is also a solution. It follows that one may solve the LSW system for large time 
iteratively by solving (|1.6p . (|1.7|) on the unit time interval < t < 1 with L(0) = 1, 
and then rescaling L{\) to 1. Hence, assuming that in the iteration one always 



DIFFUSIVE LSW EQUATION 



5 



has for some constants a, b satisfying 1 < a < b the inequality a < L(l) < b, the 
function L(t) increases linearly in t as t — ► oo. This is a quantitative statement 
concerning the phenomenon of coarsening. 

Consider now how this iteration scheme would apply to the diffusive LSW system 
(jl.7[) . (jl.lip . Denoting the solution to by c e (x,t), L £ (t), then dilation using 

the parameter A results in a function Xc £ (Xx, Xt) which is a solution to (|1.7j) . (jl.lip 
with s in (jl.lip replaced by e/X. Thus in the iteration scheme described above 
e ~ l/L(t) — ► for large time corresponding roughly to an equivalence of e ~ 
1, t — ► oo and e — ► 0, i ~ 1. 

We summarize our main results on convergence of solutions of the diffusive 
LSW system to solutions of the classical LSW system. In §4 we prove that if 
c £ (x,t), L e (t), t > 0, is the solution to (jl.7p . (jl.lip with initial data satisfying 
(jTTHl) then 

(1.13) lim / c e (x',t)dx' = / c {x',t)dx', x>0, 

Jx Jx 

limL e (i) = Lo(t), t>0, 

where co(x,t), Lo(t), t > 0, is the corresponding solution to the classical LSW 
system (jl.6[) . (jl.7l) . This is the content of Theorem 4.1. In §5 we show that 
the point-wise in time coarsening rate for the diffusive LSW system converges as 
e — > to the corresponding rate for the LSW system. Specifically, let us define 

Mf), t > o. b y 

(1.14) A e (t) = y xc e {x,t)dxj J c E {x,t)dx, t>0, 

so A e (t) is the mean cluster volume at time t > 0. Then if the initial data satisfies 
(|1.12p . the function A e (<), i > 0, is diffcrentiable for e > and 

dAM) dAJt) 
1.15) lim — -r J - = — r^, t>0. 

v ; e^o dt dt 

This is the content of Theorem 5.1. 

In [1] an upper bound on the time averaged coarsening rate for the LSW system 
is obtained by application of a rather general argument of Kohn and Otto [8] . At 
the end of §3 we extend the method of [4] to prove an upper bound on the time 
averaged coarsening rate for the diffusive LSW system (|1.7p . (jl.lip with < e < 1, 
which is uniform as e — * 0. This is the content of Theorem 3.2. 



2. Existence and Uniqueness for the Discrete Case 

We begin by giving a proof of the existence and uniqueness of solutions to the 
Beckcr-Doring system (jl.ip . (|1.2p which is somewhat different than that given in 
pp. We assume that the initial data q(0), I > 1, satisfies 

oo 

(2.i) <a(o) = 7/, i> i, £*w = p>o, 

where the 7^, £ > 1, are all non-negative. For any S > we define a space X$ by 
A5 = c\ : [0, <5] — > R : Ci(0) =71, Ci is continuous and < Ci(t) < p, < t < 8 >. 
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Now Xg is a complete metric space under the uniform norm, || • \\g. We define a 
mapping T : Xg — > Xs which has the property that the solutions to the Becker- 
Doring system give rise to fixed points of the mapping. To do this we consider 
functions v(£, t),£>l, t > 0, and define J(£, v, t) for I > 1, t > 0, by 

(2.2) J(£,v,t)=atc 1 (t)v(l,t)-b i+ i v(£+l,t), 

where ci(t), t > 0, is a given non-negative function. Let c(£,t), £ > 2, £ > 0, be the 
solution to the linear system of equations, 

(2.3) ^M=J(*-W)-J(WM>2, 
with initial and boundary conditions given by 

(2.4) c(l, t) = d(t), t > 0, c(*, 0) = 7 * , £ > 2. 

By the maximum principle there is a unique positive solution to (|2.3p . (|2.4p . The 
mapping T : Xg — > -X"$ is defined as follows: For ci G -X",? let c(£,t), < i < <5, be 
the corresponding solution of (|2 . 3[> . (|2.4[) . Then Tc\ is given by the formula 

(2.5) Tc x {t) =max|p-^&(£,t),o|, < t < 5. 

Evidently Tci £ Xg. We show that for <5 sufficiently small T is a contraction. 

Lemma 2.1. There exists S(p) > depending only on p, such that for < <5 < 5(p) 
the mapping T on Xg is a contraction. 



Proof. We have from (|2.3|) that 

oo oo ~t oo 

(2.6) J2 ic ( e '^=J2 £ ^ + ds Ahc,s)+J2 J(t,c,s). 

1=2 1=1 • J ° 1=1 

It follows that there is a constant K > 0, independent of p, S, such that 

oo oo oo 

jrX>(4t)<E*7* + *[i + p]HX^ 1/3 <#,-)IU> 

^=2 £=2 £=1 

for < t < <5. It follows from this that for 6 sufficiently small, depending only on 
p, there is the inequality, 

oo 

(2-7) ||^c(MIU<^iP, 

i=\ 

for some constant K\ independent of p and 5. 

Next we consider ci, d x <E Xg and for < A < 1, let c A .i = \c\ + (1 — A)c' x . Define 
ca(^, i), ^ > 2, < * < 6, to be the solution of ([13]), with c A ,i substituted for 
c\. Putting u\{£,t) = dc\(£,t)/d\ it follows that u\ satisfies the equation, 
(2.8) 

^ = J(^-l,«A,*)-^»«A,*)+[ci(t)-^(t)] [0«_iC A (€-l,t)-a < CA(€,t)] ) £ > 2, « > 0, 

with the initial and boundary conditions, 

(2.9) u A (l,t) = ci(t) - cUi), * > 0, u A (*,0) = 0, ^>2. 
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Let v x (£,t),£ > 2, < t < 5, be the solution to the equation, 
(2.10) 
dv x (£,t) 



dt 



J{e-l,v x ,t)-J(£,v x ,t) + \\c 1 ~c' 1 \\ 5 [a e _ 1 c x (e-l,t)+a i cx(£,t)},e>2,0< t < 6, 



with the initial and boundary conditions, 

(2.11) v x (l,t) = \\a-4Ws, t >0, u A (*,0) = 0, l>2. 

By the maximum principle v\(£,t) is positive for £ > 2,i > 0, and there is the 
inequality 

(2.12) \u x (£,t)\<v x (£,t), £>2, 0<t<6. 
From (|2~T0|) , (|2~TTj) we have that 

OO „t OO 

JT,* 1 ' 3 v x (£,t) = / d S {2 1 /3j(l )UA)S ) + ^^ + i)l/3_^/3jj ( ^ UA)S) | 



=2 



„t oo 

|C! -c'ilU / ds[2 1 / 3 a 1 c x (l,.s)+Y, [(^+l) 1/3 +^ 1/3 ]a, CA (^, S )} 



Using (|2.7[) and applying Gronwall's inequality in the previous equation, we see 
that there is a constant K 2 independent of p and S such that 



(2.13) 



provided S is taken sufficiently small, depending only on p. 
To see that T is a contraction for small <5 observe that 

oo 

\T Cl (t)-Tc[(t)\ < \ Y,t[a(Z,t)-co(£,t) 

i=i 

Hence on using (|2.6p we have that 

\Tc x (t)-Td x {t)\ < K 3 [l+p] f dsfM^-c'^l+Vy/ 3 ! f ux (£,s)d\\\, 

Jo k e=2 Jo ' 



0<t<5, 



for a constant K 3 independent of p, 5. The result follows from the last inequality, 
(l2~T2l . (|2~T3l) . □ 

Theorem 2.1. TTie Becker-Ddring system \1.2]) with initial condition satis- 

fying 112.1]) has a unique global in time solution. 

Proof. It is evident that a solution of the Beckcr-Doring system (|1.1|) , (|1 ,2[) with 
initial condition (|2.1|) gives rise to a fixed point of the mapping T of (|2.5|) . Hence 
by Lemma 2.1 the solution is unique. Correspondingly, suppose C\ is a fixed point 
of T. Similarly to (|2.6[) we have that 

00 r i f* +T ' r 

J2^,t)-c(£,t + r,)\ =- ds[j(l,c,s) 

i=i 

J(Z, c,s)\ > -K p I ci{s)ds, 

„_n ' Jt 



1=2 



8 



JOSEPH G. CONLON 



for 0<?7<(5/2, < t < 5/2, where K > is a constant independent of 5, p. Since 
Tci = c\ it follows from ()2.5p that ci satisfies the inequality, 



i-t+n 

c\(t + ?/) — c\{t) > —Kp I ci(s)ds. 

Hence if ci(0) > then ci(t) > for < t < 5/2. Conversely if ci(0) = there 
is a to > such that ci(t) = for t < to and c\(t) > for to < t < 5/2. Suppose 
to > 0. Then (|2T6j) implies that 



1=2 £=2 

Since Tci = c\ we conclude from (|2.5p that ci(t) > for < t < to, a contradic- 
tion. Thus to = whence ci(t) > for < t < 5/2 in all situations. It follows then 
from (|2.5p that the conservation law (|1.2[) is satisfied. Thus ci(t),c(£,t),£ > 2, are 
a solution of the Becker-Doring system (|1 .2|) . □ 

Next we establish existence and uniqueness for (|1.9p , (jl.lOp subject to the Dirich- 
let boundary condition c(l,t) = 0, t > 0. For a given non negative function 
ci(t), t > 0, let c(£,t), i>2, t > 0, be the solution to (|2.3p subject to the initial 
and boundary conditions, 

(2.14) c(l, t) = 0, t > 0, c(*, 0) = 7 ,, * > 2, 

where the satisfy (|2.ip with 70 = 0, p = 1. Normalizing p — p CI i t = 1 in (|1.10|1 . 
the function a(t),t > 0, is then to be chosen to maintain the conservation law 



r\ OO OO 

=2 £=1 



(2.15) ]T>(^t) = l. 

fc2 

Assuming we have a solution to (|2.3|) . (|2. 14[) . (|2.15|) it follows that 

_a 

at 

= a (t) 2 a,c(*, t) - b e c{e, t)-b 2 c{2,t). 

1=2 1=2 

We conclude therefore that ci(t) is given by the formula, 

{00 ^00 
oi?5^c(€,t) + 6ac(2,t)l / ^a t c(£,t). 
e=2 ) t=2 

Setting f = in (|2.16p yields then the identity, 

{00 ^00 
aiq ^ "tl + h 72 > / X! ai11 ■ 
1=2 ) 1=2 

Observe from (|2.16p that c\ satisfies the inequality, 

z s < ci(t) < 2{z s +9/2 1 / 3 ), t > 0. 
Hence for any 5 > we define a space X$ by 

X$ = |ci : [0, <5] — > R : ci(0) is given by (|2.17p . c\ is continuous and 
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z s < ci(t) < 2{z s + q/2 1/3 ), < t < S\. 

As previously we equip X$ with the uniform norm || • \\g and define a mapping T 
on Xg by 

(2.18) Tc l {t)=z s + \a l q^c{l,t)+b 2 c%t)\ /^a t c(£,t), 0<t<6, 



1=2 



1=2 



where c(£,t), t > 0, is the solution to 
conditions (|2.14|) . 



subject to the initial and boundary 



Lemma 2.2. There exists Sq > such that for < 5 < Sq the mapping T on X$ 
is a contraction. 

Proof. We proceed as in Lemma 2.1. We consider c\, c\ <G X$ and for < A < 1 
we put ca.i = Aci + (1 — X)c' l . Define c x (£, t), £>2 7 < t < S, to be the solution 
of ([23]), ([2TT4]) with ca : i substituted for c\. Putting u\(£,t) = dc\(£,t) / dX we see 
from (j2~Tg)) that 

/ ,1 rf°° lf°° 1 

(2.19) T Cl (t)-Tc[(t) = ^ dX r^a t c x {£,t)\ I a l9 ^ u x {£, t) + b 2 u x {2 7 1) \ 



a iq J2cx(e,t) + b 2 c x (2,t) I lj2atu X (£,t) 



1=2 



.£=2 



^aic\(£,t) 



(=2 



Evidently u x (£,t), £>2, satisfies (|2.8[) with initial and boundary conditions given 

by 

(2.20) u x (l, t) = 0, t>0, u x (£, 0) = 0, i > 2. 

Letting v x (£, t) be the solution to (|2.10[) with zero initial and boundary conditions 
corresponding to (j2~20|) it follows that (|2~T2"]) holds. 

We wish to bound the RHS of (|2.19| in terms of ||ci — c' x || 5. To do this we first 
observe that from (|2.10p we have the identity, 

^v x (£,t) = -b 2 dsv x (2,s)+2\\c 1 -c' 1 \\ s ^atc x {£,s), 



1=2 

whence we conclude that 



(2.21) 



Oil* < 25\\c 1 -c' 1 \\ s \\J2a e c x (£,-)\\ s . 



Now from 



we have that 



00 00 00 

^2a£C X (£,t) =^2aae + / eLs{a 2 J(l,c, s) +^ \a i+1 - a t J(£, c, s)| 



Defining now g(t), t > 0, by 



9 



(t)=J2a e c x (£,t), 



1=2 



it follows that there is a constant K > such that 

\g(t) - 5 (0)| <K [ g(s)ds, t>0. 
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Hence h(t) = \g(t) — g(0)\ satisfies the inequality, 



h(t) < Kt g{0) + K h(s)ds. 

Jo 

We conclude that h satisfies the inequality, 

h(t) < [e Kt - 1] 5 (0), t > 0. 

We may therefore choose <5o small enough so that g(Q)/2 < g(t) < 3g(0)/2, < t < 

8 . Now from (|2~T2"]) . (f2~2Tj) we see that 

(2.22) 



nq J2ux(£,t) + b 2 u x (2 : t)\ /J2 

aec x (£,t) 



1=2 



<Jfi%i- cilia, 0<t<5, 



for some constant K\ independent of i5o, provided < 8 < 8q. 
Next we observe from (12.81) that 



°° ft 00 

^aiux(£,t) = / ds^a 2 J(l,u\,s) + ^2^a e+1 - ae J(£,u\,s)^ 

£=2 ^° 1=2 

+ / ds a(s) ~ cj(s) (at+i - ae)aec\{£, s). 
Jo 1 e=2 

It follows then from (|2.12|) that there are constants if 2 , if 3 > independent of So 
such that 



1=2 



t=2 



1=2 



0<t<5. 



Hence, on using (|2.21j) . we see that there is a constant if 4 such that 



(2.23) 



J2a e u x (£,t)\ I J2a t c x {£,t) < KJ\\ Cl - c'^g, 0<t<8. 



1=2 



It is clear now from (|2.19|) . (|2.22p . (|2.23|) that on taking 5o sufficiently small the 
mapping T is a contraction for < 8 < Sq. □ 

Theorem 2.2. The system t2.3\) , \2.1J$ , S2.15)) has a unique solution global in 
time. 

Proof. From Lemma 2.2 the solution exists for time < t < 8. Since the only 
constraint on the initial condition is (|2.15|) at t = and on c\(i) equation (|2.17|) 
at t = we may propagate the solution beyond time 8 since these conditions also 
hold at t = 8. □ 

It is evident from (|2.16|) that for a solution to ([2~3]) . (|2.14|) . (|2.15|) the corre- 
sponding function c\(t) satisfies ci(i) > z s , < t < 00. We shall show that the 
hminf of c±(t) as t — > 00 is z s . 

Lemma 2.3. Let c(£,t), £ > 2, be a solution to i2.3\) . {2.1$ , \2.15\l . and define 
9it) by 



g(t)=J2c(e,t). 



1=2 



Then g is monotonic decreasing and lim hIX) j(t) = 0. 
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Proof. From (|2.3|) we see that g'(t) = —b2c(2, t) < whence g is strictly monotone 
decreasing. From (|2. 16|) there is the inequality. 



d(t)>z s + q/ (X 1/:i ) t , 



where (-) t is the probability measure on the integers £ > 2 with weights c(£,t). 
Using the Jensen inequality we have that 

(* 1/3 ) t < (X) 1 /" = l/g(t)^, 

whence we conclude that 

ci(i) > z. + qg(t)V a . 

Suppose now that g(t) does not converge to as t — ► oo. Then there exists 7, T > 
such that <?(i) > 7 > for t > T. We shall show that this implies 



(2.24) 



lim > c(£, t) = 00, 



1=1 



contradicting (|2.15p . 

To prove (|2.24j) we first observe that for any £ > 2 one has that 



(2.25) 



c(£,s)ds < 00, £ > 2. 



For £ = 2 the inequality (|2.25p evidently follows from the relation g'{t) = —62 c(2, t) 
and the fact that g is always positive. For £ = 3 it follows from (|2.25|) for ^ = 2 and 
the fact that g(t) — c(2,t) is bounded for all t, and we can proceed by induction to 
establish it for all £ > 2. Now using (|2.6|) we have that for t > T, 



00 „t 00 

= 1+ ds{-b 2 c{2,s) + ^[aeci(s)-bi]c(£,s)} 



> K(t) + 



ds Q7 1 / 3 g>(s), 



where by (|2.25p K(t) is bounded below uniformly in t as t — > 00. Hence (|2.24|) 
holds. □ 

Lemma 2.4. There is the limit, liminf c\{t) = z s . 

t — >oo 

Proof. We assume that liminf C\(t) > z s and obtain a contradiction. Thus we are 

t — *oo 

assuming that there exists T, 7 > such that c\ it) > z s + 7 for t >T. Let / be a 
positive function on the set {£ £ Z : £ > 2}. Then from (|2.3|) we have that 



dt 



OO OO 

£ t) = /(2) J(l, c, t) + [f{i + 1) - /(*)] ^(A c, t) 

£=2 £=2 

= I E + x ) - /(*)] a ^ ci w - - ^ ^ x )] b < I c ^ *)• 



where we put /(l) = 0. We define now /(^) by the recurrence, 

'f(£ + 1) - /(£)! 0/(2!. + 7) - \f(£) - f(i - 1)1 b e = 0, /(l) = 0, 
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for I = 2,3, £(7), where £(7) is an integer depending only on 7. For t > £(7) 
we define f(£) by the formula, 

f(i) = f(L(rf)) [^/L( 7 )] 1/3 . 

It is evident that £(7) can be chosen in such a way that there is the inequality, 



f(£ + 1) - f(£)\ a e (z s + 7) - [/(*) - /(I - 1) 
We conclude therefore that 



bp > 0, £ > 2. 



£/(*)c(*,i) > T), t>T, 



«=2 



«=2 



whence it follows that 



liminf YV /3 c(i,t) > 0. 



f=2 

Now from (|2.16|) and Lemma 2.3 we have that lim ci(t) = z s , a contradiction. □ 

t — >oo 

Remark 1. The analogue of Lemma 2.4 for the LSW equation is proved in [14] . 
Proposition 3.2. 

3. Existence and Uniqueness for the Continuous Case 

In this section we prove a global existence theorem for solutions to (jl.lip subject 
to the constraint (|1.7j) . with non-negative measurable initial data satisfying (|1.12p 
and Dirichlet boundary condition c(0, t) = 0, t > 0. We shall also assume without 
loss of generality that in (|1.7p one has p— p crlt = 1, and in p. lip that e = 1. Observe 
that the parameter L(i) in (|l.lip ought to be given in terms of the solution c[x, t) 
by the formula, 



(3.1) 



L(t) 



1/3 



x 1/ ' i c(x 1 t)dx 



c{x, t)dx 



in order that the conservation law (|1.7p holds. Thus L(t) is proportional to the 
average cluster radius at time t. 

Let Ct denote the differential operator on the RHS of p. lip and C,\ its formal 
adjoint. Thus C\ is given by the formula 

(3.2) c* t ={l + x)^ 



1 - 



L(t) 



d_ 

dx 



One can estimate the solution of (jl.lip by solving the equation dw/dt = —£*w 
backwards in time. Hence if for some T > 0, w(x, t), x > 0, t < T, is the solution 
with Dirichlet boundary condition w(0, t) = 0, we have that 



(3.3) 



w(x, T)c(x, T)dx 



w{x, 0)c(x, 0)dx. 



It is well known [7] that solutions to parabolic equations can be written as ex- 
pectation values. In particular let X(s), s > t, be the solution to the stochastic 
differential equation, 

X( S )' 1/3 1 



(3.4) 



dX(s) 



L(s) 



ds + V2 (1 + X(s)) 1/6 dW(s), 
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with the initial condition X(t) 
at 0, then 



If T Xi t is the first hitting time for the process 



(3.5) w(x, t)=E [w (X(T)) ; Tx , t > T] , x > 0, t < T, 

is the solution to the equation dw/dt = —C*w, t <T, with terminal data w(x, T) = 
wq(x) and Dirichlet boundary condition w(0,t) = 0. Wc shall show that we can 
control the solution to this terminal-boundary value problem by perturbation theory 
uniformly for x > if L(t) has a positive lower bound and for T — t sufficiently 
small. 

Lemma 3.1. Suppose that the parameter L(t), t > 0, is continuous and satisfies 
L(t) > Lq > 0, t > 0, where L < 1. Let w(x,t) be the function \3. 5\) with wq = 1. 
Then there are positive universal constants C\ , Ci such that 



Cix/VT-t < w(x,t) < C 2 x/Vf~t, < x < VT-t, d < w(x,t) < 1, x > VT-t, 

1/2 

provided t lies in the region < T — t < 6 L , where S > is universal. 

Proof. We first show the lower bound on w. Since the terminal data wo is positive 
and monotonic increasing it follows that w(x,t) > u(x,t) where u(x, t) satisfies the 
equation, 

,„ „s du . d 2 u du 

(3.6) 



du 1/c , d 2 u 

8t + V + x)/ W 



Ox 



= 0, t < T, 



with Dirichlet boundary condition u(0, t) = and terminal data u(x, T) = wq(x), x > 

0. We can estimate the solution of (|3.6[) using perturbation theory provided T— t « 

1. Let G(x,t) be the fundamental solution of the heat equation, 

(3.7) G(x,t) = 



'Aid. 



exp 



4f 



t > 0, — oo < x < oo. 



For x, y > 0, t < T, we define the kernel Kt(x, y, t) by 

(3.8) K T (x, yJ)=G (x - y, (1 + y) 1/3 (T -t))-G (x + y, (1 

We have that 



y) 1/3 (T-t)) 



dK 7 



+ (l+x) 1 / 3 



d 2 K T dK T 



dt dx 2 dx 

[a (x - y, (1 + y) 1/3 (T - *)) + [(1 + *) 1/3 - (1 + y) 1/3 ] b (x - y, (1 + y) 1 / 3 ^ - t)) } 

(z-y, (l + y)V3( T _ t) )_| a ^ + 1 , ) ( 1 + y)l/3 (T _ t )) + [ (1 + a .)l/3_ (1+y) 

6 (x + y, (1 + y) 1/3 (T - *)) }G (.x + y, (1 + y) 1 / 3 ^ - t)) , 



G* 



1/3 



where the functions a, 6 are given by 

a(x,t) = x/2t, b(x,t) = x 2 /At 2 - 1/2*. 
Consider now the function ?;(x, t), x > 0, t <T, defined by 



(3.9) 



i) 



K T {x,y,t)w {y)dy. 



Then one sees from the previous equation that 



(3.10) 



& v q _|_ n 1/3 

9* 9x 2 9a; 



<7(a;, t), x > 0, t <T, 
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where the function g satisfies the inequality, 

(3.11) \g(x,t)\<C/VT^t, 0<T-i<l, 

for some universal constant C > 0. We also have that v(0,t) = and v(x,T) = 
wq(x). Hence v{x,t) is given by the formula, 



v(x, t) = u(x, t) — E 



J g(X{s),s)ds; T x , t >T 



We conclude therefore that w(x, t) > v(x, t) — 2Cy/T — t. From (|3.9[) we can obtain 

a lower bound on v(x,t). To see this first observe from (|3.8p that 

(3.12) 

K T {x, y, t) > m G (x -y, (1 + y) 1/3 (T - i)) , \x-y\ < x/2, x > y/T-t, T-t < 1, 

for some universal rfr > 0. Since (1 + x) 1 / 3 (T- 1) < 2x 2 for x > y/T-t, T-t<l, 
we also have that 

(3.13) / G(x-y, (l + y) 1/3 (T-t)) dy > m , x > VT~ 1, T - t < 1, 

J\x-y\<x/2 V J 

for some universal 772 > 0. Hence from (|3 . 1 2|) . (|3.13|) there exists a universal 773 > 
such that v(x, t) > 773, x > y/T - t, T - t < 1. We conclude that w(x, t) > C\ > 
provided T — t < 6 < 1 for some suitable universal constant 5 > if x satisfies 
X > y/T^t. 

To conclude the proof of the lower bound we are left to consider the region 
< x < y/T-t. We write the solution of (pUP]) with zero Dirichlet and terminal 
conditions v(0, t) = 0, v(x, T) = 0, as a series 



(3.14) 



v(x, t) 

V n (x, t) 



v n (x,t), 



n=0 



ds 



dy K s (x,y,t)g n (y,s), 



{ d id 2 d 1 

.90 = .9, <?«+! = 9n - <yg t + (1 + x) 1/3 - g^j v n , n > 0. 

The recurrence relation for g n can more simply be written as 

f T f°° ( d d 2 d 1 

g n+1 (x,t) = J ds J dy j— + (1 + x) 1/3 — - — j K s (x,y,t)g n (y,s). 

It follows now in exactly the same way that we derived (|3.11[) that there is a 
universal constant C > such that 

g n (x,t) < C n (T-t) n ' 2 - 1 / 2 , n>0. 

Hence the series p,14p converges for T — t sufficiently small. Observe also that the 
sum of the <?„, n = 0, 00, is bounded by C/y/T — t for some universal constant 
C. Using the fact that 



Kt(x, y, t)dy < C min 



1, x y/T-t 
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we can bound the function v in (|3.14p . 
t + (T- t)/2. Then one has that 



Suppose x = y/so — t where t < s < 



/I roo rs . 

ds J dy K s (x,y,t) / VT^s < CJ ds j VT~ 



ds 
Cx 



ds 



/ y/{ S -t){T- S ) 



< 2C 



Vr^t- y/r = 



*0 



Cxy/2 

VT~t 



(T+t)/2 



ds 



t 



ds 



Cxy/2 
VT - t J(r+t)/2 VT - 



X, 



for some constant C . The lower bound follows easily now by obtaining a lower 
bound on the function (|3.9|) and using the previous inequality. 

We consider next the upper bound. Since it is clear that w(x, t) < 1 we need 
only restrict ourselves to x in the region < x < \JT — t. Let u(x, t) satisfy the 
equation, 



1/3 



- 1 



du 
dx 



= 0, t < T, 



with Dirichlct boundary condition u(0, t) = and terminal data u(x, T) = wq(x), x > 
0. Since the data wq is positive and monotonic increasing we have that w(x,t) < 
u(x, t). To obtain an upper bound on u(x, t) we first note that u(x, t) = P{T x ,t > T) 
where r x , t is the first hitting time at for the process (|3.4|) started at x at time t 
with L(s) = Lq. For < x < L^ 4 let r* t be the first exit time from the interval 
[0,Ll /4 ]. Evidently one has 

P(r x . t >T) <P(t* >T)+u(x), 



where u(x) is the probability that the process started at x exits the interval [0, L 
through the boundary L^ 4 . Now u(x) satisfies the boundary value problem, 



1/4, 



(1 



a/3 



d 2 u 
dx 2 



1/3 



- 1 



^ = 0, < x < Ll /4 , u(0) = 0, u{Ll /4 ) = 1. 



The problem is explicitly solvable with solution, 

u(x) 









r L o 




/ exp 

10 


— / h{z')dz 
. Jo 


dz j 


/ exp 
Jo 


— / h{z')dz 
. Jo 



dz, < x < L 



1/4 



where h(z) = [(z/Lq) 1 / 3 — l] /[l + z] 1 / 3 . Since L < 1 we conclude that there is 

a universal constant C > such that u(x) < Cx/ Lq . Hence to obtain the upper 
bound on w(x,t) for < x < yT — t it is sufficient to obtain an upper bound on 
P« t > T). 

We can do this by perturbation theory just as we did for the lower bound. 
Thus setting u(x,t) = P{r* t > T) it is clear that u(x,t) satisfies p,15|) for < 

x < l]J 4 , t < T, with terminal condition u(x, T) = 1 and boundary conditions 
u(0,t) = u(L^ 4 , t) = . Let Gjj{x,y,t) be the Green's function for the heat 
equation on the interval [0, L^ 4 ] with Dirichlet boundary conditions. Thus Go is 
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given by the method of images as an infinite series, 

oo 

(3.16) G D (x, y, t)=J2 (-l) p(m) G(x - y m , t), 

where yo = y and y m , m > 1, are the multiple reflections of y in the boundaries 
0, Lq , with p(m) being the parity of the reflection, p(0) = 0. For < x, y < Ljj 
and t < T wc define the kernel Kj>(x, y, t) by 

(3.17) K T (x, y, t) = G D (x, y, (1 + y)^ 3 (T - tfj . 
We have then that 



(3.18) 



dKT +(l + x)^ d2KT 



dt 



dx 2 



X \V3 

T ) 



dKj 
dx 



oo 

(-l) p{m) G ((x - y m ), (1 + y)^(T tj) {a{x-y m , (l+y)^(T-t)) [l- (-^) 
(1 + x) 1 / 3 - (1 + y) 1 ' 3 ] b(x - y m , (1 + y) l '\T - t))}, 



l/3i 



m=0 



just as in 



Consider now the function v(x,t), < x < L]J 4 defined by 



(3.19) 

Then one has that 

dv 
~dt ' 



v(x, t) 



K T (x,y,t)dy, t < T. 



(1 



a/3 



d 2 v 
dx 2 



1/3 



- 1 



dv 
dx 



where the function g satisfies the inequality, 

\ g (x,t)\ <c Il\' 4 Vr^t. 

Arguing as previously we thus obtain the upper bound on w(x,t) for < x < 
^JT — t provided T — t < 5L^ 2 for suitably small S independent of Lq. □ 

Lemma 3.2. Suppose that the parameter L(t), t > 0, is continuous and satisfies 
L(t) > Lq > 0, t > 0, where Lq < 1. Let w(x,t) be the function US. 5\) with 
wq(x) = a; 1 / 3 , x > 0. Then there are positive universal constants C\,C2 such that 

C lX /(T-t) 1/3 < w(x,t) < C 2 x I (T-t) 1/3 , < x < VT~t, 

dx 1/3 < w(x,t) < C 2 x 1/3 , x > VT~^t, 

1 /2 

provided t lies in the region < T — t < 5L , where S > is universal. 

Proof. The proof of the lower bound is similar to the proof of the lower bound 
in Lemma 3.1. For the proof of the upper bound however we need to make a 



modification since wq(x) 



,1/3 



is an unbounded function. Our starting point is 



as before that we wish to find an upper bound on the solution to (|3.15p . Then if 
X(s) is the stochastic process started at x at time t which is associated with the 
PDE (|3~T5l) . we have that 

(3.20) u(x,t) = E[wo{X{T))-T x . t >T] 
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< E[w (X(T));r* >t > T] + u(x) sup u(lI'\ A, 

t<t'<T v ' 

using the notation of Lemma 3.1. The first term on the right can be bounded above 
using perturbation theory as before. We shall therefore be finished if we can show 

that u(Lq /4 , A < CL\I 12 for t<t'<T, provided T — t < Li' 2 . We can use Ito 

calculus to prove this. Thus we have that 

ds + V2 ( 1 + X{s)j dW{s). 



dX(s) = 



X(s) 



Lr 



- 1 



It follows that if X(t) = x, then for t' > t, 



X(t') 



'\2 



X(s) 



L 



1/3 



- 1 

t' 



ds 



+ 2 J (l + X(s)) V3 ds + 2^2 J X(s)(l + X(s)y /6 dW(s) 



Setting t = t x t wc conclude that 



E[X(t'/\r) 2 ] = x 2 +E 



X(s) 



X(s) 



1/3 



- 1 



ds 



('At 



1 + X(s)) 



1/3 



ds 



We take now x — h\j 4 and restrict t' by t' — t < z . It follows that there is the 
inequality, 



1/2 



E[X(t' A r) 2 ] < CiLq 72 + C 2 L 1/2 £ 



-1/2 



/f'AT 
X(s) 2 ds 



for some universal constants C\, C 2 . Gronwall's inequality therefore yields, 
E[X{t' 1\t) 2 ] < CiLo /2 exp \c 2 {t' - t) / L l J 2 . 

We can now apply the Chebyshev inequality to conclude that u{L]J 4 , i') < CL^ 12 . 

□ 

Lemma 3.3. Letw(x,t) be the function (3. 5\) withwo = 1. ThenO < dw(x,t)/dx < 
C/VT^t, provided < (T - t) < SL 1 ^ 2 . 

Proof. We have already observed that < duu(x,t)/dx so we need an upper bound. 
We first prove the upper bound for < x < L 1 Q /4 /2. To do this we observe 
that w(x,t), < x < Lq^ 4 , t < T, satisfies the diffusion equation with termi- 
nal conditions wq(x) = 1, < x < L^ 4 , and boundary conditions w(0,t) = 

0, w(L^ 4 ,t) < 1, t < T. We write w(x,t) — Wi(x,t) + w 2 (x 1 t) where both wi 
and w 2 satisfy the diffusion equation. The function w\{x,t) has terminal data 
Wi(x,T) = wq(x), < x < L^ 4 , and boundary data wi(0,t) = wALq^ ,t) = 0. 
The function w 2 (x, t) has terminal data w 2 (x, T) = 0, < x < h\j 4 , and boundary 
data w 2 (0,t) = 0, w 2 (Ll /4 ,t) = w(Ll /4 ,t). 

We first consider the function wi(x,t) which we construct by perturbation ex- 
pansion. The first term in the series is the function v(x,t) of (|3.19[) . Setting 
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g(x, t) = dv/dt + Cf V, then 

oo 

wi{x,t) = v(x,t) - v n {x,t), 



n=0 



v n (x,t) 



ds 



.1/4 



dy K s (x,y,t)g n (y,s), 
d ) 

— + C t >v n , n>0. 



90 = g, 9n+l = 9n 

Just as before we see that the functions g n satisfy the inequality, 
(3.21) \g n (x, t)\ < C n (T - t) n ' 2 - x ' 2 I 4 n+1)/4 , n > 0. 

It easily follows that \dwi/dx\ < C\j \JT — t for some constant C\ provided < 
T — t < 8Ll' 2 . 

Next we consider the function W2(x,t). This can be represented in terms of a 
Green's function G(x, y,t,T),0<x,y< L^ 4 , t <T, which is the Dirichlet Green's 
function for the operator L* t on the interval [0, Lg ] . Thus if 



u(x, t) 



'1/4 



G(x,y,t 1 T)w {y)dy, t<T, 



then 



+ £*u = 0, < x < Ll /4 , t < T, 

u(0, t) = m(Lo /4 , t) = 0, u(x, T) = w {x), < x < I^ 4 . 
The function W2(x,t), t < T, has the representation, 



(3.22) w 2 (x,t) 



(io /4 ,s) ^ ( 



dG 



-1/4 



1/4 



1/:! 



To estimate W2(x,t) we compute the Green's function by perturbation series ex- 
pansion as we did before. The first term in the expansion for G(x,y,t,T) is the 
function Kx(x,y,t) of (|3.17p . If we replace G by this in (|3.22[) it is easy to see 
that \dw2/dx\ < C(T — t) / L^ 4 , < x < Lq 4 /2. The complete expansion for 
G(x,y,t,T) is given by 

OO 

(3.23) G(x, y, t, T) = K T (x,y, t) - ^v n! T(x,y,t), 



71=0 



V n ,T(x,y,t) 

9o,T(x,y,t) 

9n+l.T 

g n +i,T(x,y,t) 



ds / dy'K s (x,y' ,t)g n , T (y' ,y, s), 



g T (x,y,t) 



-+£* 

at Ux 



Kr(x,y,t), 



9n.T 



ds 



d 1 

— +C t > v ni T, n > 0, 

- 1/4 

J ° (9 1 

dy ' I dt + C * >x f Ks ^ x 'y'^^n,T{y',y,s). 
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In ()3.23|) we have written CI = C* t x to denote that the operator L\ acts on the x 
variable. Analogously to (|3.2ip we have the estimate 



(3.24) \g ntT {x,y,t)\< 



C n(r r _ t jn/2-l/2 



L 



(n+l)/4 



Glx-yJl + L, 



-1/4 



1/3 



(T-t) 



for some universal constant C. Hence the series (|3.23[) converges for < T — t < 

1 /2 

SL . We need now to differentiate the series term by term with respect to the y 
variable. Consider first the function vq,t whose derivative is formally given by the 
expression 



(3.25) 



0,T 



dy 



(x,y,t) 



ds 



-1/4 



Analogously to (|3.24|) one has the estimate, 



dgr(x,y,t) 



dy 



< 



C 



(T-t)L L 



dy' K s (x,y ,t)-^-(y ,y,s). 



Yjl G (x — y, (l + 4 /4 ) 1/3 (T-t) 



This estimate gives a nonintegrable singularity in (|3.25[) . so we need to integrate 
by parts in (|3.25|) in the s,y' variables. First we write 



(3.26) 



dy 



(x,y,t) = - 



AT+t)/2 ^ Q gT 

j ds j dy K s (x,y ,t) — (y ,y, 



s) 



ds 



-1/4 



'(T+t)/2 

We have now from (|3~T7)) . (j3~18)) that 



{ §- s + G,* 



dK 7 
dy 



(3.27) 



d 



c 



dK 7 



dy 

(f + t/') 1/3 dK T 



(y',y,s) 



L(s) 



3(1 + y) 4 / 3 dt 



V 



1 + y 



1/3 



1/3' 



d 2 K T 
dy'dy 



d 2 K T 
dsdy 



We substitute the RHS of (j3~2T)) into the second integral on the RHS of (|3~2"6"]) . 
We then integrate by parts w.r. to y' for the first term on the RHS of (|3.27p , and 
w.r. to s for the second two terms. Note that since K s {x, 0, t) = K s {x, L^ 4 , t) = 
there are no boundary terms in the integration by parts w.r. to y' . Once this is 
accomplished we can estimate dvo^/dy since we have removed all non integrable 
singularities. Thus we obtain the bound. 



(3.28) 



dv ; 



dy 



■(x,y,t) 



C 



(x-y,( 



l + L 



i/4 



1/3 



(T-t)). 



Similarly we also have that 



(3.29) 



dy 



(x,y,t) 



< 



C 



Ll /2 VT~t 



Glx-yJl + L 1 ^ 



(T-t)). 



Once we have the estimate of (|3.29p it follows by the same method as was used to 
derive ([3~2"4]) that 



(3.30) 



dg n ,i 
dy 



(x,y,t) 



< 



c n (T -ty/ 2 - 1 

r(n+l)/4 
L 



G 



(x-y, (l + Ll /4 ) 1/3 (T-t)), n>l 
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Hence from (|3.23[) we also have that 
(3.31) 



d 2 v„ 



dxdy 



(x,y,t) 



< 



C n (T - 1)"/ 2 ^ 1 /2 



-(n+l)/4 



G 



l + L 



i/4 



1/3, 



{T-t)), n>l. 



Now, just as we derived (|3 . 28[) we can see that (13.311) also holds for n 
conclude therefore that there is a universal constant C such that 



0. We 



*°(z L 4 ' 4 t s) 



< 



C 



0<x< Ll /4 /2, 0<s-t<SLl 



1/4 



1/4 



< x < 



It follows then from (|3.22|) that one has \du)2(x,t)/dx\ < C / L Q 
T-t<6L\'\ 

We consider next the case x > L,]J 4 /2. First note that we have shown that 

dw(x,t)/dx < C/Ll /A , x = L 1 /4 /2 1 T-t < 5Lj /4 . Thus if v(x,t) = dw{x,t)/dx, x > 

Lq^ 4 /2, t < T ', then it is clear that v is the solution to the terminal- boundary value 
problem, 



3(1 + x) 2 / 3 \L(t) 



1/3 



- 1 



dv 
dx 



m v = 0, t <T, x > Ll /4 /2 



3{s 2 L(i)} 1 /3 



v(Ll /4 /2,t) =dw/dx(Ll /4 /2, t), t<T; v(x,T)=0, x>L 4 ^/2. 
Let X(s) denote the diffusion corresponding to the equation (|3.32[) and for 
X(t) = x > L^ 4 /2, let T Xit be the first hitting time at L^ 4 /2. Then there is 
the inequality, 

f fTAT. t 



v(x, t) < 



C 
rV4 



E 



ds 



cxp ■ 



In view of the fact that T-t < 8Ll /2 and X(s) > L L Q /4: /2, t < s < T A r X)t , 



3[X(s) 2 L(s)] 1 /3 



V4, 



X(t) 



WC 



conclude that v(x, i) < C\jLj for some universal constant C\ if x > L ' J 2. □ 

Lemma 3.4. Let w(x,t) be the function h3. 5\) with wq{x) = x x ^ 3 ,x > 0. Then 
< dw(x,t)/dx < C/(T -i) 1/3 , provided < (T - t) < SL 4 /2 . 

Proof. We proceed exactly as in Lemma 3.3 but this time using the results of 
Lemma 3.2. Thus wc first show that dw(x,t)/dx < C / (T - t) 1/3 for x < l^ 4 j2 

and dw(x,t)/dx < C/L^ 6 for x = L^ 4 /2. For x > h\j 4 j2 wc consider the 
terminal-boundary value problem (|3.32[) but now with the terminal data given by 
v(x,T) = l/3x 2 / 3 . Since v(x,T) < Ci/L 1 / 6 ^ > L^ 4 J2, one concludes just as in 
Lemma 3.3 that dw(x,t)/dx < C/L^ e , x > L\j 4 /2. □ 

Theorem 3.1. Suppose c(x, 0), x > 0, is a non negative function which satisfies, 

poo 

< / (1 + x)c(x, 0)dx < oo. 
Jo 
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Then there is a unique solution to ll.ll\) subject to the constraint, 



xc(x, t)dx 



xc(x, 0)dx. 



Proof. Just as in §2 we define a mapping on functions ci(T) = l/ifT) 1 / 3 where 
L(T) is defined by (|3.ip . In particular, ci(0) is determined by the initial data, 
where ci(0) = l/L 1 ^ 3 for some L\ > 0. Suppose now we are given ci(T), < T < 

1 /3 

To, with ci(0) = ■ Then we solve p. Ill) with the corresponding function 

L(T), < T < T , and define a new function Ac x {T), < T < T , by the RHS of 
(|3.1|) . We define a space X of functions ci : [0,T ] — > (0,oo) which are continuous 
and satisfy ||ci||oo < l/^c/ 3 ^ or somc > 0. We shall show that the mapping A 

1/2 

leaves X invariant provided Lq < L\ is sufficiently small and To = SL for some 
universal 5, < 5 < 1. To sec this we write 



(3.33) Aci(T) 



c(x, T)dx 



x 1/3 c(x,T)dx 



Wi (x, 0)c(x, 0)dx 



W2{x, 0)c(x, 0)dx, 



where w\{x, t), t < T, is given by (j3 . 5|) with wq = 1, and W2(x,i), t < T, is given 
by (|3.5p with u>o(x) = x 1 / 3 , a; > 0. In view of Lemma 3.1 and 3.2 we have then 
that 



(3.34) 



Aci(T) < 



C 2 /"^ f 00 
— = / a; c(x, 0)dx + / c(x,0)cfa 



Vt 



a; c(x,0)dx + Ci / x 1/3 c(x, 0)<ir 

1/2 

provided T < <5T , for some universal 5 > 0. Observe now that 



x 1/3 c(a;, 0)cte < 



rl/3 ,Li/8 
^1 



1 

< - 

- 2 



c(x, 0)<ix 
x 1/3 c(x,0)dx - 



ii/8 



x 1/3 c(x, 0)dx 



Li/S 



x 1/3 c(x,0)dx, 



1 Z" 00 

x 1/3 c(x, 0)dx > - / x 1/3 c(x,0)da 

i 2 ./n 



whence we conclude that 
(3.35) 

It follows then from (|3.34|) that 

Aci(T) < 2max[C 2 ,l]/CiT 1 /3 , T < Tj/64. 
One also has from (|3.34|) that 

(3.36) Aci (T) < [C 2 /Ci + 1/CiJ/T 1 / 6 , < T < <5i /2 . 

It is clear now that there is a universal constant r\ > such that if we choose 
Lo = ?;min[Ti, 1] then A leaves X invariant. Next we show that for To sufficiently 
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small the mapping A is a contraction. To see this we proceed as in §2. Thus if 
ci ,(/]£! then 



A Cl (T) Ac[ (T) 



J\\± A[\ C1 + (1 - xy^T) 



dX\ I x 1/3 c x (x,T)dx 



dX 



(x, T)dx 



/*oo o . / r poo 

c x {x,T)dx J x 1 / 3 -^-{x,T)dxj / [y xV 3 c x (x,T)da 

where c A is the solution to (TITTTj) with l/L^) 1 / 3 = Aci(t) + (1 - A)c' 1 (t). Hence 
= dc\/d\ is the solution to the initial-boundary value problem, 

^ = C t u x - [ Cl (i) - c'M^x 1 / 3 ^}, t > 0, 

u x (x,0) = 0, x > 0, u A (0,t) = 0, i>0. 



It follows that 

/>oo 

(3.37) / u\(x,T)dx 



[ci(t) - dy(t)]dt [ dx^-(x,t) x y 3 cx(x,t), 



x 1/3 u x {x,T)dx = / \ci(t) -c[(t)]dt / Gb^(x,i)x 1/3 c A (x,i), 

where the functions u>i,u>2 arc as in (|3.33j) . Hence from Lemma 3.3 there is a 

universal constant C such that 

(3.38) 

ux{x,T)dx J j x 1 ^ 3 c x {x,T)dx 



< c 7 (T)Vr|| Cl - cilU, o < r < slI /2 , 

where y{T) is defined by 

7(T) = sup / x 1/3 c A (x,t)dx/ / x 1/3 c A (x, T)dx. 

0<t<TJo I Jo 

If we use the second equation in (|3.37[) and (|3.36|) we also conclude there is a 
universal constant C such that 



(3.39) 



cx(x,T)dx / x 1 ^ 3 ux(x, T)dx j 



x 1/3 c x (x,T)dx 



<C7(r)VrHci-</ 1 || 00 , o<r<^L, 



1/2 



Using (|3.35|) and Lemma 3.2 we see that there is a universal constant C such that 



7(T) < Cmax 



T/Lf 



1/3 



< T < 



1/2 



It follows now from (|3.38|) , (|3.39p that if we take To = v min [Lq^ 2 , L^ 5 ~\ for some 
universal constant v > then the mapping A is a contraction. Hence we obtain a 
unique solution to the diffusive LSW problem up to time To. Observe also that for 
sufficiently small L\ one has that Aci(To) - l/L 2/15 < 1/l\ /3 , whence one obtains 
global existence in time. □ 
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We show here that the Kohn-Otto argument [8] may be applied to prove time 
averaged coarsening for the diffusive LSW model (|1.7|) . (jl.l 1|) . Letting c £ (x,t) be 
the solution to the diffusive LSW system (|1.7[1 . (|f .f f p . we denote by E £ {t) the 
energy 



(3.40) 



Ee{t) 



x 2/3 c £ {x,t)dx 



From the conservation law (jTTTJ) , the quantity 1/E £ {T) 3 is a measure of the average 
cluster volume at time T. 

Theorem 3.2. Let p — p clit in \1. 7[ ) be normalized to 1, and Cq(x), x > 0, be a 
non-negative function which satisfies 



(3.41) 



(1 + x 4/3 )c (x)dx = M < 



Suppose further that c £ (x,t) is the solution of the diffusive LSW system \1. 7[ ), hl.ll]) 
with initial data c £ (x,0) = c (x), x > 0. TTien /or < e < 1, t/iere are universal 
constants K\, K 2 such that 

-3/2 



(3.42) 

provided T > K 2 M^. 



E e {tfdt 



Proof. We first show that I? e (i) is decreasing. In fact we have from (|1.11|) and (|3.1 
that 



(3.43) 



3dE 
2~dt 



e 
3 Jo 



x^^il + x/eY^c^x^dx- / x-^ s c £ (x,t)dx 



c £ (x, t)dx 



x 1/3 c £ (x,t)dx < 0. 



Next define a length scale M £ {t) by 
(3.44) M £ (t)= I x 4/3 c £ (x,t)dx, 



so that M £ (0) < Mq. It follows again from (|l.llj) that 
(3.45) 

*° x-^il+x/e) 1 / 3 ^, t)dx+ (°° x 1 / 



4 dt 3 J 



L £ {t) 



1/3 



- 1 



c £ {x, t)dx, 



where L £ (t) is given by (|3.ip with c e (-,t) in place of c(-,t). Now just as in [1] we 
have from the Schwarz inequality, 



(3.46) 



,1/3 



1/3 



2 



c £ {x, t)dx 



Le{t) 

i 

x _1 ^ 3 c e (a;, t)dx 
where we have used (|1.7[) . (|3.1I) . 



< 



xc £ {x 7 t)dx j x 1 ' 3 



c e (x, i)d:c 



x 1 ^ 3 c £ (x, t)dx, 



1/3 



2 



- 1 



c £ (x, i)da; 
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We show that for < e < 1 there is a constant C > such that 
(3.47) 



• / x~ 2/3 (l + x/e) 1/3 c £ (x,t)dx 
o 



<Ce x- i/3 (l + x/e) 1/3 c £ (x,t)dx. 
Jo 



To see this observe that the LHS is bounded above by 



(i 



2 

< 2 
+ 2 

< 2 4 / 3 e 2 

< 2 4 / 3 



i I x- 2/3 (l + x/e) 1/3 c £ {x,t)dx 
o 

i 

x- 4/3 {l + x/e) 2/3 c e (x,t)dx 



x- 7/3 {l + x/s) 2/3 c £ (x,t)dx 



x-^ 3 (l + x/e)^ 3 c e (x,t)dx 



x- 4/3 {l + x/e) 1/3 c e (x,t)dx 



x-^ 3 (l + x/e) 1/3 c £ (x,t)dx 
c e (x, t)dx 

oo 

x c £ {x, t)dx 

(1 + x)c £ {x, t)dx 



(1 + x)c £ (x, 0)dx 



We conclude then from ([3~43]) - ([07l) that 

(3.48) \dM £ /dt\ 2 < K 3 \dE £ /dt\ , 

for some universal constant K3, provided < e < 1. It is easy to see from the 
Schwarz inequality that 



(3.49) E £ (t)M £ (t) > 1, t > 0. 

The result follows from (f3T48|) . ([3^9]) and Lemma 2 of [19]. 



□ 



4. The zero diffusion limit 



Our goal in this section will be to establish the limits (|1.13p under the assumption 
(|1 . 12(1 on the initial data for the diffusive LSW system (fl~7)) . First let us 

recall how the method of characteristics can be applied to give a representation for 
the solution of (|1.6|) . For any x,t > let x(s), < s < t, be the solution of the 
terminal value problem, 



(4.1) 



dx 



L(s) 



1/3 



x(t) = x. 



We set F(x,t) = x(0), whence F(-,t) is a mapping on the non-negative reals with 
derivative 



(4.2) 



dF(x,t) 
dx 



exp 



ds 



1 



3 7 {x{sYL{s)y/\ 



The solution Cq(x, t), x, t > 0, of (|1.6p is given in terms of the initial data co(x), x > 
0, by the formula, 



(4.3) 



co(x',t)dx'= / co(x')dx', x,t>0. 

JF(x,t) 
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Lemma 4.1. Let L(-) be a continuous function on the interval [0, T] satisfying 
infL(-) > and F(x,t), x > 0, < t < T, be the corresponding mapping de- 
rived from \4- !]) ■ Suppose Co(-) is a nonnegative function on (0, oo) satisfying 
Jo°° co(x)dx < oo, and c £ (x,t), L £ (t) is the solution of U.ll\) with initial data 
co(-) and Dirichlet boundary condition c £ (0,t) = 0. Then if the functions L £ {-) 
converge uniformly to the function L(-) on the interval [0,T] as e — ► 0, one has 



lira 

£^0 



c e (x', T)dx' 



Co(x')dx', x > 0. 



lF{x,T) 

Proof. As in §3 we use the solution to the adjoint problem. Thus let C* t £ be the 
operator, 



(4.4) 



C* t e = £ (1 + x/e 



1/3 



dx 2 



1 



1/3 



9x' 



and u) e satisfy dw e /dt = — >Q £ u> 6 , t <T, w £ (x,T) = wq(x), w £ (0,t) = 0. Then 



(4.5) 



wo(x)c e (x, T)dx 



w £ (x, 0)co(x)dx. 



For any xq > we take ^(x) = 1 5 x > xo, wq(x) = 0, x < xq. Evidently 
w £ (x, 0) < 1, x > 0. Hence to prove the result we need to show that 



(4.6) 



limw £ (x, 0) = 1, x > F(.to, T), 

e— >0 



= 0, x<F(x Q ,T). 

Recall now that w £ (x,0) = P(X £ (T) > x ; T xfi > T\X £ (0) = x), where X £ {s) is 
the solution to the stochastic equation, 



(4.7) 



dX £ (s) 



Le(s) 



1 V3 1 ( \V6 



and is the first hitting time at the boundary for the process started at x at 
time t. We first show that 



(4.8) 



lim P(t x0 > T\XJ0) =x)=0, x< F(0,T). 



By our assumptions there exists £o > and Lq > such that for e < £o, inf L £ (-) > 
Lq. Hence if Y £ (s) is the solution to the equation, 

1/31 W \Va 

ds + v^efl + y e (s)/£j <W(s), 



(4.9) 



dY £ (s) = - 



1 - 



then there is the inequality, 



(4.10) 



P^t > T|X e (t) =x)< P(r x , t > T\Y £ {t) 



= x), t<T. 



Now let u e (x) be the probability that the process Y £ started at x < 5 exits the 
interval [0, S] through the boundary 0. Then one has 



u £ (x) 



where 
(4.11) 



exp 



h e (z) 



h £ {z')dz' 
(z/L ) 



dz 



exp 



h £ (z')dz' 



dz, 



1/3 



/£[1 + Z/£] 1/3 
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Choosing 5 = Lq/2, it is easy to see that on any interval [0,6'} with S' < S the 
functions u £ converge uniformly to 1 as e — ► 0. Next for < x < 6 let v £ (x) be 
defined by v £ (x) = E[t x [Y £ (0) = x\ , where t x is the first exit time from the interval 
[0,5]. Then v £ satisfies the equation, 



(4.12) 



1 /3 

(l + x/ej v"{x) + 1 - (x/L ) 1/3 v' s (x)=l, < x < 6, 



with Dirichlet boundary condition v s (0) = v £ (S) = 0. It follows from the maximum 
principle that 

v £ (x) < x I [l - (S/L ) 1/3 ] , < x < S. 
We conclude therefore that 

(4.13) P(Y s {s) hits before time t \Y £ (0) = x 



> u £ (x)-P(r x >t\Y £ (0) = x^ 

> u £ {x)-x/t 1 - {S/L ) 1/3 . 

Note that if x << t then the RHS of (|4. 1 3[) is positive as e — * 0. Now on setting 
e = in (|4.12p we see that the time for the classical system to exit the interval 
[0, S] is T x , where 



dz/\l- (z/L ) 1/3 



Hence we should have that 

(4.14) lim P (t x >T x + V \Y £ (0) = x) = 0, 



for any r\ > 0. In order to prove (|4.14[) we first show that lim e ^ v e( x ) = T x . From 
(|4.12p we have that 



v' e {x) = A £ exp 



h £ (z)dz 



exp [— J x h £ (z')dz'] dz 
e(l + z/ey/ 3 



where A £ is given by the formula 



A £ = 



i: h £ {z')dz' 



exp 

\<z<x'<8 e(l + 2/e) 1 / 3 



dzdx' 





r' 




/ exp 


— / h £ (z)dz 


dx 


Jo 


Jo 





Letting g(z) — 1 — (z/Lq) 1 ^ 3 we have that 



* e*v[- £ h e {z')dz'} 



:(l + z/ey/ 3 



dz 



exp 



1 



1 



exp 



g(x) g(o) 

Similarly we have that 



h £ (z')dz' 



x h £ (z) 

o g{z) 

x Sfiz) 



h £ (z')dz' 



dz 



■ exp 



h £ (z')dz' 



dz. 



1 



5(0) 



L 



0<z<x'<5 



9^z) 
9i z )' 



■ exp 



I s dx ' 1 1 






Jo g&) 1 J 


exp 










- [ h £ (z')dz' 


dzdx 


11 


J z 







h £ (z')dz' 



exp 



dx' 



h £ (z')dz' 



dx'. 
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We conclude that 



g{z) 2 

6 

exp 





■ , a 


f 5 dx' 




r' 


exp 


— / h £ (z)dz 


Jo 9(x') i 


1 J exp 


— / h £ (z)dz 




. JO 




Jo 



dx' 



exp 



h £ {z')dz' 



{h{x-z) 









- [ h £ (z')dz' 


dx' 1 J exp 


- f h £ (z')dz' 


Jo 




Jo 



dx 



:'}dz, 



where H{y) is the Heaviside function, H(y) = 1, y > 0, H(y) = 0, y < 0. It follows 
easily from the previous expression that lim e ^o v' e {%) = ^/9( x ) uniformly in any 
interval [0, 6'] with 5' < S. Hence lim s ^o v e (x) = T x for < x < S. 

We can in a similar way estimate w E (x) = B[tJ|5^(0) = x] . In fact w E satisfies 
the equation, 



(4.15) 



z(l + x/e) 1/3 w"(x) + [l - (x/L Q ) 1/3 } w' £ (x) = 2v £ (x), < x < 5, 



with Dirichlct boundary condition w s (0) = w £ (S) = 0. Proceeding as in the previ- 
ous paragraph we see that lim e _ > o w' £ (x) = 2T x /g(x) uniformly in any interval [0, 5'] 
with 5' < S. It follows that lim e ^o w s (x) = T x for < x < S. Finally we conclude 
from the Chebyshev inequality that <|4.14[) holds. 

From (|4.10|) and (|4.14|) we see that for x satisfying < x < L a /2 then (|4~8l) 
holds provided T > T x . We show now that it holds for all x < F(0, T). To do this 
let f x , t be the first hitting time for the process X e (s) with X £ (t) = x > Lq/4 to hit 
Lq/4. Putting f = T x fi we see as in Lemma 3.2 that 



E 



XJt' At) 



2E 



Xe(s) 



2eE 



Le(s) 

(l + X £ (s)/s) 1/3 d S 



1/3 



ds 



It follows that there is a constant C{Lq) depending only on L$ such that 

ft' Af 



E 



x £ (t'A?y 



< x 2 + C(L ) E 



X £ (s A f) 2 ds 



< t' < T. 



GronwaU's inequality therefore yields, 



(4.16) 



E 



X £ (t'Af) 2 < x 2 exp [C{L )t'} , 0<i'<T. 



Next, applying Ito's lemma to [X £ (t') — X (t')] we see from f|4. 7[) that 



[x £ (t')-x (t')y 



[X £ (s) - X (s)} 



1/3 



X (s) 
L (s) 



1/3 



ds 



2V2e 



[X £ (s) - X (s)] [1 + X £ (s)/e} 1/6 dW(s) +2e [ [1 + X £ (s)/e] 

Jo 



1/3 



ds . 



With r as in (|4.16p we then have that 
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(4.17) 

E{{X E (t' A?)-X (t' A?)} 2 } <C(L ) 



sup \L £ (s) - L (s)\ + e 

0<s<T 



2/3 



: exp [C(L )f] 



/" r ~i 

+ / E{[X e (s At) - X {s At)} 2 dsj , 0<t'<TAT„ 



for a constant C(Lq) depending only on Lq, where T x is defined by Xq(T x ) = 0. 

Again from the Ito lemma we have that 

(4-18) 

x.m-*m = [ {xM}" 3 -{TMr d ' + ^Sl [1 + w^wi 

Now by the Kolmogorov inequality we have that for 77 > 0, 

{'Af \ q„ T /-TAf 

[l + X e (s)/e} 1/e dW(s) 



P 



sup 




\0<t'<T 


Jo 



[l + X e (s)/s} 1/3 ds 



We also have that 



P sup 

\ 0<t'<TAT x 



t'Af 



LJs) 



1/3 




for a constant C(Lq) depending only on Lq. It follows then from (|4. 1 6[) . (|4.17p . 
P~Tg)l that for any 77 > 0, 



(4.19) 



lim P I sup 

\0<t'<TAT x Af 



X e (l/) - Xotf) 



> v = 0. 



We finish the proof of ([18]). For x < F(0,T) note that < T. We take ry in 
(|4~19|) so that < L /4 and consider the event {\X £ (T X A r) - X (T X A r)| < ry}. 
We evidently must have f < T x for this event. Since then |io/4 — Xo(r)| < ?/ the 
stopping time f must be close to the time required for the classical trajectory to go 
from x to Lq/4. Now by starting the diffusion X e at time t = f and using (|4.10p . 
(|4~T4|) we conclude that ([478]) holds. 

To complete the proof of the lemma we need to show (|4.6[) . For x > F(0, T) one 
has T x > T whence the event { | X e (T x A f ) — X (T x A f) | < 77} has high probability. 
If f > T then r 2 , i0 > T so we consider the situation f < T. Since \L / A — X$(t)\ < ?] 
in this case, the stopping time f must again be close to the time required for the 
classical trajectory to go from x to Lq/A. Using our previous argument we can 
estimate the time taken for the diffusion X e started at time t — f at Lq/A to hit 
0. For e small this is close to the time for the classical trajectory, whence r Kj o > T 
with high probability. We conclude that 



limP^o >T\X e {0>) 



)=1, x>F(0,T). 

This completes the proof of (|4. 6|) for xq = 0. The argument for xq > is similar. □ 



Lemma 4.2. Let c e {x,t), L e {t), t > 0, be a solution of the system \1. 7\ ), hi. 11}) 
with initial data co(x) satisfying cq(x) > 0, x > 0, (1 + x)co(x)dx < 00, 



1 f°° 1 

xco(x)dx > — / xco(x)dx = — , 
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J °° xco(x)dx = 1. TTiera for any eo, T > £/ie set of functions L e (-),0 < e < £o, 
on t/ie interval [0, T] is an equicontinuous family. 

Proof. Let us first consider the classical LSW system (|1.6p . (|1.7[) and let A(i), i > 0, 
be the mean particle volume as in (11.141) . Then one has by Jensen's inequality that 
L(t) < A(i) and it is also easy to see that A(t) is an increasing function. For x > 
let Wq(x) be defined by 

/>oo 

wo(x) = / co(x')dx' . 

J X 

From (|4.3p we have that A(i) = 1/^0(^(0, i)) and one can further see that F(0, t) < 
t. Noting that A(0) is given by 

/>oo / />oo 

A(0) = / xco(x)dx I j co(x)dx, 

we see that 

/>oo 

/ xco(x)dx > 

Ja(o)/2 2 Jo 2 

whence u>o(A(0)/2) > 0. Hence A(-) is bounded above in the interval [0, A(0)/4], 
and a-fortiori L(-) is bounded above in [0,A(0)/4]. Next we obtain a lower bound 
on L(-) in this interval. To do this we use the formula (|3.ip for L(t) which, if we 
denote by w(x, t) the LHS of (|4.3|) . is the same as 

i r°o 

L(t) 1/3 = - / x~ 2/3 w(x, t)dx /w(0, t). 

Since F(x, t) < x + t we have that 

1 f°° 1 f°° 

- / x^ 2/3, w(x,t)dx > - x~ 2/3 w (x + t)dx 

3 Jo 3 J 

> a 1/3 w (a + <), a > 0, t > 0. 

We have then that 

£(f) 1/3 > A(i)a 1/3 w (a + t), t > 0, a > 0. 

If in this last inequality we set a = A(0)/4 we see that L(-) is bounded below on the 
interval [0, A(0)/4] by a positive constant. This argument can be easily extended, 
using the fact that A(-) is an increasing function, to conclude that A(-) is bounded 
above on any interval [0, T] and £(•) is bounded below by a positive constant on 
the interval. 

The boundedness below of L(-) on [0, T] implies the continuity of the functions 
A(-), £(■). To see this note that for 6 > 0, 

< Aity 1 - A(t + oT 1 < wo (F(Q, t)) - w (F(S, t)) . 

In view of (JOJ) we have that < F(S,t) - F(0,t) < S. The continuity of A(-) 
follows now from the continuity of wo(-). To show continuity of L(-) we note that 
since 

/>oo 

/ w (F (x,t))dx = 1, 
Jo 

there exists for any r\ > an x n > such that 

/■OO 

(4.20) / x- 2/3 w Q (F(x, t)) dx < rj. 
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We consider the integral 



(4.21) 



where 



x' 2/3 [w (F(x, t)) - wo (F(x, t + 5))} dx 



2.1 ., 



x~ 2 / 3 [wo (F(x, t)) - wo (F(g(x, t, S),t))} dx, 



\g(x,t,6)-x\ <6\l + (x/L )V 3 



Using the continuity of wo(-) and (|4.2p we conclude that 5 can be chosen sufficiently 
small that the integral in (|4.21[) is bounded in absolute value by r\. Hence from (|4.20|) 
we have that 



x~ 2/3 [w {F(x, t)) - w {F(x, t + 5))) dx 



< 377. 



This last inequality and the continuity of A(-) implies the continuity of -£*(•). 

Next we wish to apply the previous argument to the system (|1.7[) . Let 
A e (t) be as in (| 1 . 14[) , whence it follows that A e (t) is an increasing function of t and 
A e (0) = A(0). From ([4~5]) we have that 



.22) 



/>oo 

A e (T) _1 = / w £ (x, 0)c (x)dx, 
Jo 



where w £ (x,T) = 1 and w £ is a solution to the adjoint equation as in Lemma 4.1. 
It is easy to see from (|4.7[) that 

w e (x, 0)>P {r xfi > T\X e (0) = x) , 

where X £ (s) is the solution to the stochastic equation, 

(4.23) dX £ (s) = -ds + Vte (1 + X £ (s)/s) 1/e dW(s). 

We shall show that there is a positive constant "f(T,Eo) depending only on T and 
£0 such that 

(4.24) P(T Xi0 >T\X e (0) = x)>i(T,e ), x > 2T, 

provided < £ < Eo- To see this we argue as in Lemma 4.1. Thus from Ito's lemma 
applied to (|4.23[) we have that 



E[X £ {thT xfi ) 2 } =x 2 -2E 



tAT m n 



X £ (s)ds 



+2eE 



l + X £ ( s )/e 



1/3 



ds 



On using the inequality, 

(1 + z) 1 / 3 < 1 + z 1 / 3 < 1 + e- 1 / 3 + e 5 / 3 z 2 , 

we see then that 



E [X £ (t A r xfi ) 2 } <x 2 + 2t(e + e 2/3 ) + 2e 2/3 £; 
It follows from the Gronwall inequality that 
(4.25) E [X £ (t A T K , ) 2 ] < \x 2 + 2t(e + e 2/3 ) 

We also have from (|4.23[) that 



X £ (tAT X:0 ) -x + t 



X £ (s) 2 ds 
2e 2 / 3 t 
[l + X £ ( s )/e] 1/6 dW( S ), 



exp 
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whence one has that 
i] 2 p( sup |-X,(f At-,,0) 



0<t<T 





fTAT x , 




1/3 


) < 2eE 




l+X s {s)/e 


ds 




Jo 







Using now (|4.25|) and the Holder inequality we conclude that 
(4.26) 

lup \X E (t A 7-3,0) - x + t\ > T) ) < 4Te 2/3 \e 2 + \x 2 + 



2T(e 



- e 2/3 )l exp 



2 e 2 / 3 j 1 



Choosing 77 = x/2 in (|4.26[) we see that there is an 17 > 2T such that the RHS of 
(|4.2Gp divided by r] 2 does not exceed 1/2 provided x > Xx, < e < Sq. Furthermore 
there exists Et > such that we may take x^ = 2T if £ is in the region < e < Et- 
We conclude that 

(4.27) P(t x , > T\X e (0) = x) > 1/2, x > 2T, < e < e T , 

> 1/2, x > x T , < e < e - 

To complete the proof of (|4.24[) then we need to show that 

(4.28) P {t x . > T\X e {0) =x)> 7(T, e ), 2T < x < x T , e T < e < e . 

This can be demonstrated using perturbation theory. Choose an integer N and 
points xj = 2JT/NJ = 0, 1, ...,N + 1. Now define p h j = 1, AT by 



inf (Probability X e {t) with X E (0) = Xj 

exits the interval [ay— 1, ay+i] in time greater than T/N} 



Then the LHS of (|4.28[) is bounded below by pi ■ ■ ■ pn ■ To estimate the pj we choose 
TV large enough so that the Dirichlct Green's function for the interval [ 
can be expanded in a converging perturbation series as in Lemma 3.1. The integer 
TV can be chosen dependent only on T because er < £ < £o- One can thus show 
that the solution of the diffusion equation corresponding to (|4.23|) with initial data 
1 and Dirichlct boundary conditions is bounded below by a constant depending 
only on T at x = Xj, t = T/N. Hence we obtain a lower bound on pj and thus 
have established (|4.28j) . 

We have now from P~2"2"j) . P~2"4]) that 



(4.29) 



A £ (T) < l/7(T,e ) / c (x)dx, 

I JIT 



whence A e (T) is uniformly bounded above for < e < £q provided T < A(0)/4. 
Arguing as before we can obtain a uniform lower bound on L e (t), < t < A(0)/4. 
To do this we write 

/•oo 2 /*oo roc 

/ x 1/3 c e {x,t)dx = - I x~ 2/3 / c s (x',t)dx', 
Jo 3 J Q J x 



whence for any a > 0, 



L £ {t) 1/3 > A £ (t)a^ 3 / c e (x',t)dx' 
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Now in the same way as we obtained (|4.29|) we see from this that there is a positive 
constant if (A(0)) depending only on A(0) such that 

/>oo 

(4.30) L £ (t) 1/3 > K(A(0)) / c {x)dx, 0<t< A(0)/4, < e < e Q . 

JA{Q)/2 

We have shown that the functions A e (-), L £ (-),0 < e < £o, are uniformly bounded 
above and away from zero in the interval < t < A(0)/4. Next we show that they 
are equicontinuous. To do this we use a formula analogous to ()4.22|) . Thus for <5 > 
we write 

/>oo 

(4.31) A e {T) - A E {T + 5)- 1 = / w E (x, 0)c o (x)dx, 

Jo 

where w £ (x,t) is a solution to the adjoint equation as in Lemma 4.1 with terminal 
data specified at t = T . The terminal data is given by 

w £ {x,T) = 1 - P (X E (s) > 0, T <s<T + S\X £ (T) = x) , 

where X e (s) satisfies the stochastic equation (|4.7|) . It is easy to see that 

(4.32) w £ (x, T) < 1 - P (X £ (s) > 0, < s < S\X £ (0) = x) , 

where X £ (s) satisfies the stochastic equation ()4.23|) . We can use (|4.26[) then to 
estimate the RHS of (|4~3"2"]) . Thus on setting T = S and 77 = VS + x/2 in (j4^S|) we 
see from (|4.32[) that w £ (x,T) satisfies the inequalities, 

(4.33) w e (x,T) < 1, 0<x<AV6, 

w £ (x, T) < C(e )6[l + x 1/3 ]/x 2 , x > 4v^, 

where C(eo) is a constant depending only on e - 

We shall show that if w s (x,T) satisfies (|433)l then the RHS of ([431]) is small 
for small S, uniformly in e, < e < eq. To see this we first consider the case where 
< e < (5 1 / 4 . Let X £ (s) satisfy (|4.7|) and Y £ (s) the corresponding deterministic 
equation, 

Y e (s) 1 ^ 



(4.34) 



dY £ (s) 



1 



ds, < s < T. 



We obtain an estimate on the size of the difference X £ (T) — Y £ (T) when X £ (0) = 
Y £ (0) = x, which is similar to (|4.26[) . Letting t x be the first exit time from the 
interval (0,oo) for the process X £ (s) started at x at s = we have as in the 
derivation of (|4.25p that 

ftAr x 



E [X £ (t A r^o) 2 ] < x 2 +2t(e+e 2 ^)+2t/L^+(2e 2 / 3 +2/L^)E 

where Lq is a lower bound on L £ (-). We conclude that 

(4.35) E [X £ (t A r xfl ) 2 ] < [x 2 + 2t(e + e 2 ' 3 + 1/Lj /3 )] cxp [2(e 2 ' 3 

From (g3]), ([4T3^)1 we have that 

(4.36) 



X £ (s) 2 ds 



E 



{[X £ (tAT x )-Y £ {tAT x )] 2 } 



2E 



tAr x 



[X £ ( s ) - Y e (s)] [XAs) 1 / 3 - Y £ ( S )^ 3 ] L^sfl^ds 

ftAT x 



2eE 



[l+X £ (s)/e} 1/3 ds 
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< 

ri 



provided Y e (s) > 0, < s < t. Setting y(t) to be the LHS of ([436]) we see that 

\dy/dt\ < 2y(t)/L 1 Q /3 Y £ (t) 2 / 3 + 2eE{[l + X £ (tAT x )/e} 1/3 } 

< 2y(t)/Ll /3 Y e (t)^ + 2(e + e 2 ' 3 ) + 2e^E [X £ (t A r x ) 2 } . 
Integrating this last inequality we conclude that 

(4.37) E{[X e {tl\T x )-Y e {tl\T x )f} 
^ dt'{2(e + e 2 / 3 ) + 2e 2 ' z E [X e (t' A r x f] } exp I J 2/ Ll /3 Y £ (s) 2/3 ds 

Using P~3"5)) we see from P~5T)) that X £ (t A t x ) differs from Y £ (t A t x ) by 0(e 1/3 ). 
Recall now that the function w £ (x, 0) in (|4.3ip is given by the formula, 

(4.38) v) e (x,0) =E[w £ (X s {T),T);t x > T\X e (0) = x] , 

where w £ (x,T) satisfies (|4.33[) . It is easy to see that the RHS of (|4.31[) is small 
for small 6 uniformly in e, < e < eo, if we replace X £ (T) in (|4.38|) by Y e (T). In 
fact let a £ have the property that if Y e (0) = a e then Y e (T) = 0, whence w e (x, 0) = 
0, x < a e . Now from it follows that if Y s (0) > a £ + (5 1 / 24 then Y e (T) > 6 1/24 . 
We conclude therefore that w e (x,0) satisfies the inequalities, 

w e (x, 0) = 0, x < a e ; w e (x 1 0) < 1, a £ < x < a e + S 1 ^ 24 , 

(4.39) 

U7 e (a:,0) < KS 11 / 12 , x>a £ + 5 1 / 24 , 

for a constant K depending only on eq. Since co(-) is an L 1 function (|4.39p implies 
that the RHS of (|4.31[) is small for small S, uniformly in e, < e < Eq. 

We can extend the argument in the previous paragraph to estimate the actual 
function w e (x,Q) of (|4.38[) by using (|4.37[) . In fact one sees using the Chcbyshcv 
inequality and w e (-,T) < 1 that 

(4.40) w e ( Xl 0) < K6 1/12 , x > a £ + S 1/24 . 

Here we are using the assumption e < S 1 ^ 4 . Next we show that w e (x,0) is small if 
x < a £ — 8 1 / 24 . To see this first note that for such an x then Y e {t) with Y e (0) = x 
satisfies Y e {t) = for some t = t £ satisfying t e < T — 5 1 / 24 . Now from (|4.37jl one 
has then that X £ (t e A t x ) < 0(S 1 ^ 12 ) with high probability. We need to show that 
this implies that t x < T with high probability. We consider the probability that 
the diffusion X £ (t) with X E (t e ) = x hits before time t E + S 1/24 . This is greater 
than the probability of the diffusion Y £ (s) defined by ()4.9j) with Y £ (Q) = x hitting 
before time 5 1 / 24 . Let Is be the interval [0,5 3 / 48 ], and u £ {x) be the probability 
that Y £ (s) with Y £ (0) = x exits Is through 0. As previously u £ (x) is given by the 
formula, 

^3/48 5 3 / 48 z 

u £ (x) = / exp — / h £ {z')dz' dz I exp — / h £ {z')dz' 



dz, 

where h £ (z) is defined by (|4.1ip . Next let v £ (x) be the expected time for the 

diffusion Y £ (s) with Y £ (0) = x to exit Ig. Then there is the inequality, 

(4.41) 

1 -p(y £ (s) hits before time S 1/24 \Y £ {0) = x) < 1 - u £ (x) + S^ 1/24 v £ (x), x G I s . 
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Comparing the RHS of (|4.41|) to the formulas obtained after (|4.14[) we conclude 
that if < x < 6 J / 96 then the RHS of (|44Tj) is bounded by 0(S 1/32 ). Combining 
this with the Chcbyshcv inequality applied to (|4.3T|) we conclude 

(4.42) w e {x, 0) < K6 1/48 , 0<x<a £ - S 1/24 . 

The inequalities (j440|) . (|442)) now imply that the RHS of ([431]) is small for small 
5, uniformly in e, < e < S 1 ' 4 . 

To complete the proof of the equicontinuity of the functions A £ (T) we need to 
consider the case S 1 ^ 4 < e < £q. Our goal will be to show that 

(4.43) w e (x,T-6 1/3 )<K6 1/24 , x > 0, 

for some constant K depending only on £q, Lq. It follows that w E (x,0) is also 
bounded by the RHS of (|443]) . whence the RHS of (j43"T|) is small for small 5, 
uniformly in e, 5 1 / 4 < e < £q. To establish (|4.43j) we shall use perturbation 
theory to solve the terminal-boundary value problem for w £ (x,t) on the domain 
< x < L , T - S 1/4 <t <T. The terminal data satisfies ([3733)) . The boundary 
data at x = is zero and we may estimate the boundary data at x = Lq from 
(I433)l . (14371) . Thus since t > T — S 1 / 4 we have as before from the Chebyshev 
inequality that 

(4.44) w e (L ,t) < K5 1 / 4 , T - 8 1/4 <t<T, 

again for a constant K depending only on Eq, Lq. 

We proceed in a similar way to Lemma 3.3. Thus we write w e (x, t) = wi^^x, t) + 
W2,e(x,t) where wi >e has zero boundary data and uu2 ie has zero terminal data. 
To estimate u>i. e , let Go(x,y,t) be the Dirichlet Green's function for the interval 
[0, Lq], where Go is given by the formula (|3.16|) . Similarly to (|3.17|) we define 
K £ ,T(x,y,t) by 

(4.45) K e<T (x, y, t) = G D (x, y, e(l + y/ £ ) 1/3 (T - i)) , < x, y < L , t < T. 
Consider now the function 

(4.46) v E (x,t) = [ ° K e , T (x,y,t)w E (y,T)dy, t <T. 
Then one has that 



dv e u . ,i/3 d 2 v e rr x ~~\dv e , . 



9v s , (A , , ^1/3 

"(*)■ 

where g e may be estimated as in Lemma 3.1. Observing that 

\e(1 + x/e) 1/3 -E(l + y/E) 1/3 \ < \x-y\, 
we conclude that there is a constant C depending only on e,Lq such that 

\g e (x,t)\<c/y/e(T-t), t<T. 

We conclude that 



(4.47) \w 1 . e (x,t)\<v e (x,t) + C^(T-t)/E. 

To estimate W2, e we note that W2, e has the representation, 

W2, e {x,t) = E w e (Lo,T Xtt );T Xtt < T, X £ (t x j) = L \X e (t) = x 
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where T X} t is the first exit time from [0, Lq] for the diffusion X e (s),s > t, with 
X e (t) = x. Hence from (|4.44|l we conclude that 

(4.48) w 2 , e {x,t) < KS 1 ^, T-8 1/A <t<T. 

The inequality (|4.43|) follows now from (|4.46|) . (|4.48|) since e > <5 1/4 . This completes 
the proof of the equicontinuity of the functions A £ (T), < e < Eq, in the interval 
< T < A(0)/4. 

Next we show equicontinuity of the functions L £ (-) 1 < e < £q, in the same 
interval. Analogously to (|4.20|) we show that for any 7] > there exists x v > such 
that 

/>oo />oo 

(4.49) / dx x~ 2/3 c s (x',T)dx' < 77, 



for all e, < e < Eq. Observe that for A > 0, 
(4.50) 

pOC pOC pOC pOO 

/ dx c £ (x',T)dx'= dzc (z) dx P (X E (T) > x; r 2j0 > T\X E (0) = z) . 

J A Jx JO J A 

From (|4.35p there is the inequality, 

P (X £ (T) > x; T,, > T\X e (0) = z) < min [l, K{z + l) 2 /x 2 ] , 
for some constant K depending only on Lq, £q, T. Hence 

poo 

(4.51) / dx P (X e (T) > x; r zfi > T\X e (0) = z) < K x min [z + l,(z+ l) 2 /A] 

J A 

for a constant K\ depending only on Lo,Eq,T. It follows easily from (|4. 5 1 [) that 
the RHS of (|4.50p can be made arbitrarily small by choosing A sufficiently large. 
Now the existence of x v satisfying (|4.49j) follows. 

To complete the proof of equicontinuity of L £ (-) we need to show that 

i>2x rl />oo 

(4.52) / dxx~ 2/3 [c e (x',T)-c e (x',T + 5)]dx' <rj 

JO Jx 

for sufficiently small 5 > 0, provided < e < £q. As in (|4.3ip we write 

/>oo />oo 

(4.53) / [c E (x, T) — c e (x, T + 6)] dx = / w e ^ a (x, 0)co(x)dx. 



Thus to prove (|4.52[) wc need to get estimates on the functions w e . a which are similar 
to the estimates we obtained for a = 0, but which are uniform for a satisfying 
< a < 2x v . This is a straightforward extension of the method we used for the 
case a = 0. 

We have proved equicontinuity of A e (T), L e (T), < e < £q, in the interval 
< T < A(0)/4. We wish to extend this now to arbitrary values of T. To see this 
suppose L e (-) is bounded below by Lq > in the interval [0, T]. Then there is a 
constant A depending only on Lq, T,So such that 

/>oo />oo 

(4.54) / dx c e (x',T)dx' < 1/8, < e < e . 

J A Jx 

This follows from (|4.50|) and (|4.51|) . To show that L e (-),A E (-) remain bounded 
beyond the interval [0,T] we note that 

(4.55) / dx c E (x' ,T)dx' + ^-t- / c £ (x' , T)dx' > 1/2. 



A e (T)/2 Jx * JA e (T)/2 
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From (|4.54|) . ()4.55|) we see that there is a constant K > depending only on 
Lq,T, £o and A e (T) such that 

/>oo 

/ c E {x' ,T)dx' >K, 0<e<e O - 
We proceed now as previously using the same methodology as we used to prove 

□ 

Theorem 4.1. Let c E (x,t) 1 L E (t), < x,t < oo, 6e the solution to the diffusive 
LSW problem \1. 7[ ), with initial data Cq(x) satisfying 

/>oo 

(1 + x)co(a;)(ix < oo, / x co(x)dx = 1. 
o Jo 

Denote by co(x,t), L(t), the solution of the LSW problem ll.6\) . {1. 7[ ) wii/i £/ie same 
initial data. Then there are for all x, t > the limits, 



(4.56) lim / c E {x' 7 t)dx' = / cq(x' , t)dx' , 

e ^°Jx Jx 

KmL e (t) = L(t). 

e^O 

Proof. From Lemma 4.2 we have that the functions L E (t), < e < Eq, are equicon- 
tinuous on any finite interval [0, T]. Hence there is a subsequence ej, j = 1,2, 
with linij^oo e.j = such that L e (t) converges uniformly on the interval [0,T] as 
e — > through the sequence {sj} to a continuous function L(t). By Lemma 4.1 
it follows that the first identity in (|4.56|) holds, where cq(x, t) is the solution to 
(|L6]1 . Since (fTTfj) also holds for c e (-,t) , it follows from ([43H]) that (fTT]) must hold 
for co(-,t), whence cq(x, t) is the solution to the LSW problem (|1.6p . (|1.7|l . By 
uniqueness for the solution to the LSW problem, we can then conclude that (|4.56|) 
holds as e — > through the reals. □ 

5. Convergence of Coarsening Rate 

Finally we wish to show that the rate of coarsening for the diffusive LSW problem 
(|1.7j) . (jl.lip converges as e — > to the rate of coarsening for the LSW model (|1.6|) . 
(|1.7p . To do this we will prove that 



f°° d f°° 

J Ce{x,T)dx=—J c (x,T)dx, 



(5-1) lim-^ 

e-+o dl 

where c e and Co are as in Theorem 4.1. Evidently (|1.15p follows from Theorem 4.1 
and (|5~T|) . Observe that from ([4T2} . (|4T5]| the derivative on the RHS of (|5~T|) exists 
provided the function Co(-) is continuous at x = F(0,T). We show next that the 
derivative on the LHS of (|5 . 1 1) exists for all e > 0. 

Lemma 5.1. Suppose the initial data Cq(x) for the diffusive LSW problem |1.7[ ). 
il.ll]) satisfies the conditions of Theorem J±.l. Then for e, T > the function 
J c £ (x,T)dx is differ entiable w.r. toT. 

Proof. Let w e (x, t, T),t<T,x> 0, be the solution of dw £ /dt = -£$ ;S w e , w e (x, T, T) 
1, w e (0, t, T) = 0, t< T, where L* t e is given by (|4~4]) . Then we have 

/■ oo 

c E (x,T)dx = / w E (x, 0, T)co(x)dx. 
Jo 
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We shall show by perturbation theory that for t < T such that T — t is sufficiently 
small, the function w e (x,t,T) is differentiable w.r. to T and the derivative is 
a bounded function. This will prove the result since the function v e (x,t,T) = 
—dw e (x,t,T)/dT satisfies dv e /dt = —C* te v £ , t < T, with boundary condition 
v s (0,t,T) = 0. 

Observe that 1 — w £ satisfies the diffusion equation with zero terminal data 
and boundary data 1 at x = 0. We shall show how to construct this func- 
tion using perturbation theory. We first restrict ourselves to some finite interval 
< x < e. Let wi tC (x,t, T), t < T, < x < e, be the solution of dw\ t£ /dt = 
-Cl e w h£ , wi, e (x,T,T) = 0, wi, E (0,t,T) = 1, wi te (s,t,T) = 0, t < T,0 < x < e. 
Just as in (|3.22[) u>i. e can be represented in terms of the Dirichlet Green's function 
G £ (x,y,t, s) for the interval. Thus 



(5.2) 



w ltS (x,t,T) = e 



T QQ 

ds — — (x, 0, t, s). 



dy 



As in (|3.23| we can represent G £ in a series expansion. 



(5.3) 



G s (x,y,t,T) 

V n ,e,T(x,y,t) 

go,e,T(x,y,t) 

9n+l,e,T 



k £ ,t(x, y,t)-y j V n ,e,T{x, y, t), 



ds 



9e,T(x,y,t) 



dy' K e>s (x, y', t)g n , £ , T (y', y, s), 
d 1 



9n,e,T 



d 

-Q t + A*e } «n,e,T, n > 0, 



where K £ t is given by (|4.45p . The function g n , £ ,T is given by the recursion formula, 



(5.4) g n +i, e ,T(x,y,t) 



ds J d y' \jf t + ^t,e^ K ^A X 'y'^)9n,e,T(y'y,s). 



Now as in p.24p there is a universal constant C > such that 
(5.5) 



C n (T - iW 2 _ 1/2 
\g n ,e,T&V,t)\ < .J,y 2 G{x-y,2e{T-t)), n > 0. 



£ (n+l)/2 

Hence the series (|5 .3[) converges for T — t < e/C. If we formally differentiate the 
RHS of (|5.5p w.r. to T we are led to expect the inequality, 



(5.6) \dg n ,e,T{x,y,t)/dy\< 



C n (T-t) 



i/2 - 1 



G{x-y,2e(T -t)), n>0. 



£ (n/2+l) 

It is easy to see that (|5.6[) holds for n = 0. For n = lwe use (|5.4[) . Writing 

AT+t)/2 i-T 

(5.7) gi, £iT = ds+ ds, 

Jt J(T+t)/2 

one can easily see that the derivative of the first integral on the RHS of (|5.7[) w.r. 
to y is bounded by the RHS of (|5.6p for n = 1. To estimate the derivative of the 
second integral we need to integrate by parts. We argue as in (|3.27p . Thus 

1/3 



(5.8) 



ds £ > y 



dK, 



e.T 



dy 



{y',y,s) 



v 



Le(s) 



1 



d 2 K t 



e.T 



dy' dy 
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(l + y'/e) 1 / 3 dK £ 



3e(l + y/e) 4 / 3 ds 



1 



(l + y/e)V3 



a 2 ^ 



9s 9y 



We substitute the RHS of (|5 .8(1 into the second integral in (|5.7[) as the expression 
for dgo ye: T(y',y,s)/dy. We then integrate by parts w.r. to y' for the first term, 
and s for the last two terms in (|5.8[) . It is easy to see from this that (|5.6[) holds for 
n = 1. Hence (|5.6[) holds by induction for all n > 0. It follows then from (|5.3p that 

C n (T — t) n / 2 

(5.9) |^„, £ ,t(z, J/, t)/Sy| < ^„ /2 + ; G(a; - y, 2e(T - t)), n > 0. 



Note that to obtain (|5.9[) for n = we need to use the representation (|5.8|) and 
integrate by parts. Hence if (T — t) < e/C the series in (|5.3|) for the derivative 
dG E (x,y,t,T)/dy converges, and we conclude that 

(5.10) \dG e (x, y, t, T)/dy\ < S— G(x - y, 2s(T - t)), 

V £ (J _ *J 

for some constant C. Now from (|5.2[) . (|5.10p we see that Wi £ (x, t, T) is differentiable 
w.r. to T and is given by the formula 

(5.11) ^ (x I t,T)=e^(x J 0,t i r). 

We return to consideration of the function w e (x,t,T). If ^Q(s) is the diffusion 
process associated with C* E then 

l-t« e (a:,t i r)=P(T a ., t <Tpr e (t)=a;), 

where t Ij( is the first hitting time at for X e (s),s > t. Suppose that < x < e 
and let Ti >X) t be the first exit time from the interval [0,e] for X £ (s),s > t, with 
X e (t) = x. If X s (ri tXj t) = e denote by T2, x ,t > Ti,x,t the first hitting time at e/2. 
The density dfi x j(s), s > t, associated with T2, x .t is defined by 

P (T 2 , Xtt < T ; X £ (T hx , t ) = e) = J dfi x , t (s). 
We can use this density to write 1 — w £ in terms of the function Wi t£ . Thus 

(5.12) l-w e {x,t,T)=w h£ (x,t,T) + J dfi x ,t{s) Wl , e {e/2,s,T) 

d/J, x ,t(si)diJ, £ /2, Sl (s 2 )u;i, e (e/2, s 2 , T) H . 



/t<S!<S 2 <T 

It is evident that the series converges and is term by term differentiable w.r. to T, 
provided (T - 1) « e. □ 

Let us write p £ (x,t,T) = —dw £ (x,t,T)/dT where the function w £ (x,t,T) is 
given in Lemma 5.1. Then p £ (x,t, •) is the density for the exit time r x t at of the 
diffusion X £ (s) given by P~T|) with X £ (t) = x. 

Lemma 5.2. There are positive constants c\, G\ > such that ifT — t = C\E then 
p £ satisfies the inequality. 



p e (x,t,T) < — cxp [-(a;/4e) 5/3 
£ L 



x > 0. 
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Proof. The result follows from Lemma 5.1 provided < x < e. For x > e let T x ,t 
be the first hitting time at e/2 for the diffusion X e (s) with X £ (t) = x. Then from 
the formula (|5.12|) it will be sufficient to show that 



(5.13) 



P(r x , t < T) < Cy exp -{x/Ae) 



5/3 



For any non-negative integer N let In be the interval [2 N e, 2 N+1 e], and Z £ (s) the 
diffusion process started at Z e (t) = x € In which satisfies the stochastic equation, 



(5.14) 



dZ e (s) 



-ds 



2e{\ + Z e {s)/e) 1/6 dW(s). 



We denote by T\, Xt t the exit time from In for the process Z £ (s) with Z £ (t) = x. 
Evidently one has P(r Xt t < T) < P(ri Xi t < T). We can generate the density 
for Tx tX) t by perturbation theory as in Lemma 5.1. The series converges provided 
(T - t) /2 N / 3 e « 1 If we take x ~ 2 N+1 /' 2 e then one sees that provided T-t = c x e 
for some sufficiently small universal constant c\ there is the inequality. 



P(ri, x ,t < t + ae) < Ct cxp \~{2 N /\0f' z 



where C\ is also a universal constant. The inequality (|5. 13|) follows. 



□ 



From Lemma 5.2 we see that the integral of the function p e (-,t,T) is bounded 
independent of e as e — > provided T — t ~ e. We shall show that the integral is 
in fact close to 1 for T — t = 0(1) as e — > 0. 

Lemma 5.3. There exists Sq > such that if S satisfies < 5 < Sq then one can 
find e(5) > for which the following inequality holds: 



(5.15) 



1 



p e (x, t, T)dx 



<S 1 / W , T-t = 5, 0<e<e(S). 



Proof. We first show that (|5.15[) holds in an averaged sense. Thus let < r\ < 8/2. 
We shall see that there exists e(8,rj) > such that for < e < £(5,rf) there is the 
inequality, 



(5.16) 



1 

1 



t+s+ v 



2r/ J t+S - V 



dT 



p(x, t, T)dx 



< 



s 1 / 4 . 



Let r x< t be the first hitting time at for the diffusion process X s (s) associated with 
C* e , where X s (t) = x. Then (|5. 16|) is the same as 



(5.17) 



1 - — / 



< 



Now (|5.17p will follow if we can estimate P(r Xt t > T) sufficiently accurately for 
T — t ~ 5. In particular we show that P(T x ,t > T) ~ 1 if x is slightly larger than 
T — t (depending on 5 and e) and P(r x , t > T) ~ if x is slightly less than T — t. 
To do this we proceed as in Lemma 4.1. We consider the situation when x is large. 
For Z e (s) satisfying (|5. 14|) let t x be the first exit time for the diffusion Z £ (s) with 
Z e {t) = x from the interval [x/2,2x\. Then P(r Xit < T) < P{t x < T). Now one 
has 



Z e {s A t x ) = x- [s Ar x -t} + 



2e [1 + Z e (s')/e} 1/6 dW{s'). 
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It follows that provided (T — t) < x/4 there is the inequality, 

2 



(5.18) 



P(t x <T)< 



2e[l + 2x/e} 1/3 (T-t) 



(5.19) lim P ( T Xyt > t- 



Hence one can choose e(6) so that if e < e{8) the integral in (|5.17p over x > 85 is 
negligible. For x < 86 we see just as in the derivation of (|4.14|) that 

+ x[l + KS 1/3 ]j = 0, Um P(t x j < t - rf + x) = 0, 

for any rf > 0, where if is a constant independent of e and 5. The inequality (|5.17|) 
follows from (|5~T8j) and ([57T9]) . 

To obtain the pointwise estimate (|5 . 1 5|) we use the fact that p £ (x,t,T) satisfies 
the equation dp e /dt + CI £ p £ = 0, p £ (0,t,T) = 0. Let c e (x,s),s > t, satisfy the 
adjoint equation dc s /ds = C Si£ c £ , c e (0, s) = 0, s > i, c £ (x, t) = 1, x > 0. Then one 
has that 



(5.20) 



p E {x, t, T)dx 



p e (x, s,T)c £ (x, s)dx, s > t. 



From p. lip we see that c £ (x, t) is given by the expectation value, 
(5.21) 



c £ (x, s) = E 



cxp ■ 



2 

9e 



l + X £ (s')/e) 



-5/3 



where X e (s') is the diffusion process satisfying 



(5.22) dX £ (s') = -<l- 



Xe(s') 



[Le(s') 



1/3 



111 




T x ,s < 



+ V2l (l + X e (s')/e) 1/6 dW(s'), 



s < s, 



with X £ (s) — x. Note that the process X £ (s') is running backwards in time. The 
stopping time t x , s is the first hitting time for the process on the boundary x = 0. 
From (|5.21[) we see that c £ (x, s) < 1, s > t. Hence by taking T — s ~ e we can 
conclude from Lemma 5.2 that the LHS of (|5.20p is uniformly bounded as e — > for 
any fixed t < T. Let p £ (x, t, T) be the density for the diffusion Z £ which corresponds 
to the density p £ (x,t,T) for X £ . Using the time translation invariance of Z £ we 
have from the inequality (|5. 16[> for p' £ that 



(5.23) 



1 

1 

2-q 



T-8+2ri 



ds 



T-S 



p' £ (x, s, T)dx 



< 



for < e < e(S, rj). Now from equation (|5.20p for p' £ one has that 



(5.24) 



p £ (x,T - S,T)dx = — 



T-S+2r) 



ds 



T-S 



p' £ {x, s, T)c' £ (x, s)dx, 



where c' £ is given by an expectation similar to (|5.2ip . Using the fact that the integral 
on the RHS of (|5.24[) is concentrated at x ~ 5 and (|5.2ip for c' e , we conclude from 
(|5.23p that one can choose e(S) such that for < e < e(S) one has the pointwise in 
time estimate 

/>oo 

|1- / p' e (x,T - S,T)dx\ < S 1 ^ 4 . 
Jo 
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The inequality (|5.15|) follows if we can show that p' £ (- ,T — 6, T) and p £ (-,T — 5, T) 
are close. 

We first compare p e (x, t, T) and p' £ (x, i, T) when T — t = c\e as in Lemma 5.2. To 
do this we write p e (x, t, T) as a sum over walks on the numbers 0, 1/2, 1, 2, 2 N , 
where 2 N e ~ 1. For a walk X(n), n = 0, 1, 2, ... we take X(0) > 0. If < X(n) < 
2 N then X(n + 1) can be cither of the neighbors of X(n). If X(n) ~ 2 N then 
X(n + 1) = 0. Finally let r be defined by X(t) = and X(n) > 0, n < r. Then 
for x = 2"e for some integer n with — 1 < n < N we have that 



(5.25) p e (x,t,T)= Yl 

{walks X(-):X(0)=x/e} 



S(s - t)6(s T - T) 



SQ<S1<S2..<St 



ds ds 1 ...ds T ] [ p s {eX(n-l),eX{n),s n -\,s n ) 



Here p E (x,x' ,s,s') is the density at x' at time s' for the diffusion X £ (-) of (|4.7p 
with X e (s) = a; exiting the interval I x where I x = [x/2, 2x] if x = 2 n e with < n < 
N, I x = [0, e] if x = e/2. For x = 2 N e then x' = and p E (a;, x', s, s') = p e (x, s, s'). 
Just as we derived (|5.10p we see that there are positive universal constants G\ , Ci 
such that 



(5.26) 



p B (x,x',s,s') < — 



■ cxp 



-C 2 x 5 / 3 /e 2 / 3 ( S '- S ) 



s — s 

It follows that we may assign transition probabilities on the walks X(-) so that 



(5.27) p E (x,t,T) < — e\ exp 



5/3 



|X(0) = Sc/£ 



for some positive constants C, 77. This last inequality gives another proof of Lemma 
5.2. To compare p E and p' e we use the representation (|5.25[) and an interpolation 
formula. Thus we have 
(5.28) 

Pe(x,t,T)-p' e (x,t,T) = J dX Y 

{walks X(.):X(0)=x/e} Sq<Si . 



S(sq— t)S(s T — T)dsodsi 



r 

Y\ \Xp e (eX{n- l),el(fi),s n _i,s n j + (1 - X)p' e (eX(n- 1), sX(n), s„_i, s n 

n=l 

T 

Y, [pe(eX(n-l),eX(n) 7 /)^£l(n- l),eX(n), 



/ [A/9 £ (£X(7j - l),£X(n), s„_i,s n ^ + (1 - X)p' £ (eX{n - l),eX(n), s n -!, s n 
Now using the perturbation method of Lemma 5.1 we see that 



C1X 1 / 3 

(5.29) \p e (x,x',s,s') - p' e (x,x',s,s')\ < — exp 



-C 2 x 5 / 3 / £ 2 / 3 ( S '- S ) 



provided x = 2 n e for any integer n, —l<n<N. Since 2 N e ~ 1 the inequality 
([5T29]) also holds for x = 2 N e. This follows from ([Q6]) . Just as we derived (f5T27|) 
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we have from (|5.29|) that 



\p e (x,t,T)-p' e (x,t,T)\ < -^i^X^^exp 

£ U=i 

whence it follows that 



5/3 



n=l 



k(0) = ar/e , 



(5.30) 



\p £ (x,t,T) - p' e (x,t,T)\ < ^exp [-C 2 (x/e) 



\5/3 



provided T — t = c\E and x < 0(1). For x > 0(1) the inequality (|5.30[) follows 
from Lemma 5.2. 

We shall now use (|5.30[) and (|5.20[) to compare the integrals of p e (-,t,T) and 
p' e (-, t, T) when T — t = 5. In fact we have 

(5.31) J °° p e (x,t,T) - p' e (x,t,T)dx = 

J °° [p s (x, T - ae, T) - p' e (x, T - cie, T)\ c s (x, T - as)dx 

+ Jo°° Pei x i T ~ c i £ i T ) [ce(x,T - cie) - c' s (x,T - cie)] dx, 

where c' e is the analogous expectation for the diffusion Z e which corresponds to c e 
for X e . Since c e (x, T — c\e) < 1 it follows from <|5.30[) that the first integral on the 
RHS of ()5.3ip is small. We shall show the second integral is also small by obtaining 
pointwise estimates on the differences c' £ (x, s) — c £ (x, s) for s <~ 5 + t and x > 0. In 
fact, in view of Lemma 5.2, we may restrict ourselves to obtaining an estimate on 
the difference when x = 0(e). 

To carry this out we consider the representation for c! e analogous to the repre- 
sentation (|5.21[) for c e . Let Y e (s') be the diffusion process satisfying 
(5.32) 



3 



ds'+VTe (l + Y E (s')/e) 1/6 dW{s'), s' < 



with Y E {s) = x. If we allow the diffusions X s (s') of ([ST2"2")) and Y E (s') of (|Q2"j) to 
be driven backwards in time by the same Brownian motion dW(s') then it is clear 
that if X £ (s) = Y e (s) = x then X e (s') < Y s (s'),s' < s. Thus if t x>s is the first 
hitting time at for X € (s') with X £ (s) — x and t' x s is similarly defined for Y e (s'), 
then t X)S > t' x s . Now similarly to (|5.2ip we have the representation, 



(5.33) 



c' e {x,s) = E 



exp 



s 2 
9e 



(l + Y e (s>)/e)- 5/3 d S i 



; r' < t 



If we compare (|5.21[) , (|5.33[) and use the fact that t x>s > t' xs , X £ (s') < Y E (s'), s' < s, 
we see that c e (x, s) < c' £ (x, s). We also have that 



c' e (x,s)-c £ (x,s) = E 
+ E 



exp 
exp 



(i + y e ( s ')/e)" 5/3 ^' 



9e 

' l(l + Y E ( s ')/e)- 5/3 ds' 



exp 



1 



—— ds' \ ; t x ,„ < t 



We conclude from the previous inequality that 
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ds' 



(5.34) 

< c' E (x,s)-c £ (x,s) < P« s < t, r KiS > t)+E ^ -^—^L——. ; Tx & 

+ E A ^ [1 + X £ ( s ')/e]" 5/3 - [1 + F £ ( S ')/ £ ]" 5/3 ds' ; r X)S < t 

We shall estimate each term on the RHS of (|5.34p . Note that for the purposes of 
doing this estimate we may replace L e (s'), t < s' < t + 6, on the RHS of (|5.22p by 
L Q = inf t < s /< t+5 L £ {s'). 

We consider the first term which is the same as P (t' x s < t) — P (r x , s < t). We 

shall estimate this by obtaining an upper bound on p(^t' x s < tj and a lower bound 

on P (t XiS < t). We have now that for < x < S 2 , 

(5.35) P (t' x s <t) < P {Y £ {s') exits [0, 5 2 ] through 5 2 \Y E (s) = x) 

+ P (Time to exit [0, 5 2 } larger than 8/2 \Y £ (s) = x) . 

Letting u e (x) be the first probability on the RHS of (|5.35|) we have that 



(5.36) e(l + x/e) 1/3 u' £ '(x) + 



l + |(l + .x/ £ )- 2 /3 



u' e (x) = 0, 
< x < S 2 , u e (0) = 0, u £ {8 2 ) = 1. 



The solution to (|5.36[) is given by the formula, 



(5.37) 



u £ (x) 



s- 



h £ (x')dx'/ / h £ (x')dx', 



h £ (x) = (1 + x/e) 



-2/3 



exp 



3 /-, , / \2/3 

-- (1 + x/e) 



Observe that u £ (x) ~ 1 for x > 0(e). We estimate the second probability on the 
RHS of (|5.35|) by calculating the expected time to exit [0,(5 2 ]. Letting t x be the 
time to exit and v £ (x) = E[t x ] then v £ (x) satisfies 



(5.38) -eil + x/e^^v'^x) 



v' £ (x) = 1, 
< x < S 2 , v £ (0) = v £ (S 2 ) = 0. 



We can write v £ {x) as an integral w.r. to the Dirichlet Green's function G £ (x,y) 
for the equation (|5.38|) . Thus we have 

r S 2 



V £ (x) 



where G £ {x,y) is given by the formula, 
(5.39) 

" r S 2 



G £ (x,y)dy, 



G £ (x,y) 

if x < y, and 
(5.40) 

G £ (x,y) 



h £ (x')dx' 



h £ (x')dx' 



r f s 2 




■ rv 


/ h £ (x')dx' 




/ h £ (x')dx' 


J X 




Jo 



I eil + y/ef^hM h £ (x')dx', 



:(l + y/e)^ 3 h £ (y) h £ {x')dx' 
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if x > y. In (|5.39|) and (|5.40p the function h £ {x) is given by (|5.38|) . Observe that 

upon integration by parts, 

(5.41) 

h £ (x')dx' = e{l+y/e) 1/3 h £ {y)-e{l+8 2 /e) 1/3 h £ {5 2 )-\ f (l+x' /e)- 2/3 h £ (x')dx' . 

3 Jy 

Hence we see from (|5.39j) that G £ (x,y) < 1 if x < y. Using the fact that (1 + 
x/e) 1 ^ 3 h e {x) is a decreasing function, we similarly see from (|5.40p that G £ (x,y) < 1 
for x > y. Thus v e (x) < 8 2 , whence by Chebyshev the second term on the RHS of 
(1 5 . 3 5 [) is smaller than 26. 

To obtain a lower bound on P(t x ^ s < t) observe first that 



(5.42) p(i 



< 



t) >P \X e {s') exits [0, VS\ through VS \X e (s) = 

x P (Time to exit [0, 2V6] is larger than S \X £ (s) = \/~8 



Letting u e (x) be the first probability on the RHS of (|5.42|) we have that 



(5.43) u £ {x)> h s (x')dx' / J h £ (x')exp |^ 



2/.'! 

'5 



1/3' 



dx . 



It is evident that the difference between the RHS of ([Q7]) and the RHS of ([5T3]) 
converges to zero as e — > provided < x < S 2 /2. To estimate the second 
probability on the RHS of (|5.42j) we let t x be the time to exit [0, 2y/5] and put 
v E (x) = E\t x \- Arguing as we did previously we have that 

(5.44) ]imv £ (VS) 



dx 



^ l-(x/L )V 3 ' 



Similarly if we set w £ (x) = E[t 2 ] then 



(5.45) 



lim ui £ (v<5) 



2V5 



dx 



^ i-{x/L y/ 3 



It follows from ([STi3]l . ([5^5)1 that the second probability on the RHS of ([5TI2]) 
converges to 1 as e — > 0. We conclude that the first term on the RHS of (|5.34[) is 
bounded by 4<5 as e — > 0, uniformly in x for < x < 5 2 /2. 

We turn to the second term on the RHS of (|5.34|) . Letting t x be the exit time 
of X s (s') with X £ (s) = x from the interval [0, 2y/S], then if we put 

ds' 



v £ (x) = E 
one has that v £ (x) satisfies 

(5.46) -e{l + x/e) 1/3 v'^{x) - 



X £ {s) = x 







1 - 


(i) 





v' £ (x) = 1/x 2 / 3 , 



< x < 2V5, v £ (0) = v £ (2y/S) = 0. 
Hence v £ {x) is given in terms of the Dirichlct Green's function for (|5.46j) by 



Ve{x) 



G £ (x,y) y 2 ' 3 dy. 
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Arguing as we did for the Green's function (|5.39|) we see that for S sufficiently small 
G £ (x,y) < 2, whence v E (x) < 125 1 / 6 . To bound the second term on the RHS of 
(|5.34[) we need to add to it paths for which t x < S and X e (r x ) = 2y/5. Such a path 
can eventually wander back to x = y/5 and then we can estimate by v £ (sfS) the 
contribution to the integral of X £ (s')~ 2 / 3 until the path again hits 2y/S. Thus the 
second term on the RHS of (|5.34() is bounded by 



(5.47) 



3L 



1/3 



v £ (x) + v £ (V5) + P(r^ < 8)v e (V5) 



.fe=2 



for some constant K. Note the last term on the RHS of (|5.46p is to take account 
of the integral of X £ (s')" 2/3 when X e (s') > VS. We are also using ()5.44j) . (|5.45|1 to 
see that P(r < 5) < 1/2 for small e. 

We consider the final term on the RHS of (|5.34[) . Let t x be the time taken for 
the diffusion X s (s'),s' < s, with X e (s) = x to exit the interval [0,2v5], and let 
t' x be the corresponding time for Y e (s'),s' < s. Since X £ (s') < Y £ (s'), s' < s, it 
follows that for a path X £ (s'), s' < s, with X £ (s — t x ) = 2vS then t' x < t x . Hence 
the final term on the RHS of (|5.34[) is bounded by 





' 2 


E 






_9e 




' 2 


+ E 






9s 


- E 


2 




9e 




' 2 


+ E 






9e 



(5.48) E ^ / [1 + X £ {s')/e) 5/3 ds' ; r x > 8,X £ (s - t x ) = 

[1 + X e ( s ')/e]" 5/3 ds' ; X £ (s - t x ) = 2VS 

[1 + Y £ (s')/er 5/3 ds' ; Y £ (s - t' x ) = 2^6 

■ J\l + X £ (s')/e]- 5/3 ds' ; t 2V - S s < t\X £ (s) = 2VS 
Letting u £ (x) be the first term in (|5.48|) we see that 



(5.49) -e{l + x/e) 1/3 u'^{x) 







1 - 









1/3 9 



— (1 + x/e)~ 5 / 3 P (t x > 5, X £ (s - t x ) = 0) , 



< x < 2VS, u £ (0) = u £ (2V5) = 0. 
Next we put v £ (x) = E[t x ; X £ (s — t x ) = 0]. Then v £ (x) satisfies 
(5.50) 



:(l+x/e) 1/3 v'^(x)- 



1/3 

j- ) 



-2/3 



v' s (x)=p(x s (s-t x )=0), 



< x < 2VS, v £ (0) = v e (2VS) = 0. 



We can see from (|5.37p that P( y X £ (s — t x ) = 0) << 1 if x > 0(e). Hence from our 
estimates on the Green's function (|5.39[) we conclude from (|5.50|) that v £ (x) < Ce 
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for < x < 2^/6, where C is a constant. Now using the Chebyshev inequality to 
bound the RHS of ([5^|) we see that u £ (x) < Ce/5, < x < 2y/5. 

Next we consider the second and third terms in (|5.48[) . The sum of these two 
terms is given by u £ (x) — v £ (x) where u £ {x) satisfies 



(5.51) -e{l + x/e) 1/3 u^(x)- 
2 







1 - 









1/3 



^(l + ^)" 2/3 



u' E {x) 



9e 



l + x/e)- 5 / 3 P(X e (s - t-x) = 2V5), 



< x < 2V6, u £ (0) = u £ (2VS) = 0, 



and v £ {x) satisfies 



(5.52) -e(l + x/e) 1/3 v',!(x) - 



1 + ^1 + ^/e)- 2 / 3 



-(l + x/e)- 6 ' 3 P (Y e (s-r> 



< x < 2V5, v £ (0) = v £ (2VS) = 0. 



It is easy to see from (|5.37[) that 

\P(X e (s - t x ) = 2V6) - P(Y £ (s - t x ) = 2^5) I < Cs 1/3 . 

We may also easily estimate the difference in the Green's functions (|5.39[) for (|5.5ip 
and (|5.52j) to conclude that \u £ (x) - v £ (x)\ < Ce 1/3 for < x < 2\fl. 

We estimate the last term in (|5.48|) in a similar way to how we estimated the 
second term on the RHS of (|5.34[) . We see that it is bounded by e 2 ! 3 /8 1 / 3 . We 
have therefore appropriately bounded the third term on the RHS of (|5.34[) . whence 
the result follows. □ 

Theorem 5.1. Let Cq(x), Co(x,t), c £ (x,t) be as in Theorem J^.l and F(x,t) be the 
mapping defined after J^. J[ ). If the function Co(-) is continuous at x = F(0, T) then 
JO) holds. 



c £ (x, T)dx 



p e (x,0,T)c (x)dx, 



Proof. We use the fact that 
d 

dT 

where p £ is as in Lemma 5.3. We first show that 

(5.53) lim / p e (x,0,T)c (x)dx = c (F(0,T)) lim / p £ (x, 0, T)dx. 
To do this we use the representation, 

(5.54) Pe(x,0,T) =E[p(X e (t),t,T) ; r xfi > t \X £ (0) = x) , 

where X £ {s) is a solution to (|4.7|) with X £ (0) = x and t x $ is the first hitting time 
at 0. We choose now t — T — c\E so we can bound p £ (-,t,T) from Lemma 5.2. It 
follows from Lemma 5.2 that for x > F(Q,T) then 

(5.55) p £ (x, 0, T) < - {exp [-(.t/8 £ ) 5 / 3 ] + P (X £ (t) < x/2 ; t Xj0 > t)} . 

It is evident that the first term on the RHS of (|5.55[) is small. We shall show that 
the second term is also small provided x >> 1 + 0{T). To do this let t x be the exit 
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time for X £ (s) with X e (0) = x from the interval [x/2,2x] and define Y £ (s) 7 s > 0, 

by 

,1/6 



t2e{l+X e (s')/e) ' dW(s'). 



Observe now that for any X,M\^ £ (s) defined by 



M\, e (s) = exp 



s > 0, 



(f + X e { S ')/s) 1/3 ds' + XY £ (s) 
is a Martingale. Hence for any a > there is the inequality, 

p( sup Y e (s)>a) < p( sup il/ A , £ (s) > exp -eA 2 (l + 2x/e) 1/3 t + \a ) 



< exp 



:\ 2 (l + 2x/e) 1/3 t- \a 



Minimizing the RHS of the last inequality w.r. to A we conclude that 
p( sup \Y £ (s)\ > aj < 2 exp -a 2 /4e{l + 2x/e) 1/3 ] 

0<s<t 

Let X £ {s) be the classical solution to (|4.7|) with X £ (s) — x, i.e. the solution 
to the deterministic equation when we eliminate the stochastic term. Then from 
Gronwall's inequality we have that 

sup \X E (s A t x ) - X c e {s A t x )\ < sup |y e (s)| exp[Ct/x 2 / 3 ], 

0<s<t 0<s<t 

for some constant C. If we use the fact that X^{s) > x — s, we may conclude from 
the last two inequalities that 



P (X £ (t) < x/2 ; r xfi >t)< P{t x <t)<2 exp -x 2 cxp(-Ci 1/3 ) / 64e(l + 2xle) x 

provided x > At. Thus we obtain an estimate on the RHS of (|5.55[) which falls off 
exponentially in x/e, and we can easily extend this to all x > F(Q,T) + 5 for any 
5 > 0. 

Next we consider the case < x < F(0,T), in which case we use the bound 
p £ (x, 0, T) < CP (r X fi > t |^e(0) = x) je for some constant C, which follows from 
(|5.54[) and Lemma 5.2. We use the comparison (|4. 10[) and let t x be the exit time 
from the interval [0,(5] for the process Y e (t) with Y e (0) = x. Then for a function 
f(x), if we define v s (x) by 



f(Y £ (t))dt\\Y £ (0) = x 



v E (x) = E exp 
then v E [x) satisfies 

(5.56) e(l + x/e^v'^x) - [l - (x/Lo) 1 / 3 ] v' e {x) + f(x)v e (x) = 0, 

< x < 5, v £ (0) = v £ (S) = 1. 

Observe now that the boundary value problem 

rju'^(x) — au'^x) + \u v (x) =0, x > 0, u n (0) = I, 
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has solution u^x) given by the formula 

a — \/ a 2 — 4?yA 



u v (x) = exp 



2r, 



x > 0, 



where we assume ij,a,X > and Ar/A < a 2 . If 77 — > then u v (x) converges to 
exp [Ax/a], corresponding to the fact that the deterministic path starting at x takes 
time x/a to hit 0. It is clear that if we choose rj = e{\+5 / e) 1 / 3, and a = 1 — (S/Lq) 1 ^ 3 
then v e (x) = u v (x) satisfies (|5.56[) with f(x) > A and i> £ (0) = 1, v e (6) > 1. We 
conclude therefore that provided ArjX < a 2 there is the inequality, 



E e 



< u n (x), < x < S, 



whence 

P(t x >t)< e' xt u v (x), t > 0. 

If we optimize now w.r. to A we conclude that P(t x > t) satisfies the inequality, 

(5.57) P(t x > t) < exp [ - (ta- x) 2 /4rjt], < x < at. 

It follows that if < x < 5 and 5 is sufficiently small, depending on T, that 
lim p e (a;,0, T) = 0, uniformly in x. 

We can apply the arguments of the previous two paragraphs to show that for any 
S > then lim p £ (x, 0, T) = uniformly for x in the interval < x < F(0, T) - 5. 

In fact we have the inequality, 



(5.58) P{r xfi > t \X e (0) =x)<P^ sup 

0<s<t-?j 



\X s (sAt x )-Xc(sAt x )\> p5/2) 

P(Ty > V), 



+ 



sup 

0<y<(l+p)5/2 



where p, 77 are small positive constants. The exit time t v in the second term on 
the RHS of (|5T55|) refers to the exit time from [0, 6} for Y £ (s) with Y s (0) = y. The 
exit time t x in the first term is the time taken for X £ (s) with X e (0) = x to exit 
the interval [<5/2,a], where a = sup{x(s) : < s < T} and x(s) is the solution to 
P~Tj) with x(T) = 0. For the solution to gH]) with x(0) = x < F(0, T) - 8, let rj 
be defined by x(T — r/) = (1 — p)S/2 for some small p > to be determined. Then 
X £ (t — ?y) < 5/2 for small e, whence (|5.58[) holds for sufficiently small e. We can 
choose now p > small enough so that 77 > (1 + 2p)S/2. In that case both terms on 
the RHS of (|5.58|) go to zero exponentially in 1/e 2 / 3 as e — > 0. We have therefore 
proved (|5.53p . 

We complete the proof of the theorem by showing that 



(5.59) 



lim 

e->0 



p £ (x, 0, T)dx = exp 



ds 



lo x(s) 2 /3i( s )l/3 

where x(s) is the solution to (|4.1|) with x(T) = 0. To do this we write as in (|5.20l) 



p e (x,0,T)dx 



p e {x, T — 5, T)c e {x, T — 8)dx. 
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Since by our previous argument p £ (x, T — 5, T) is concentrated at x ~ S, (|5.59|) will 
follow from Lemma 5.3 if we can show that 



where x(s) is the solution to (|4.1|) with x(T — 5) = x ~ S. This follows as before 
using the representations (|5.21| . (|5.22| for c £ (x, s) and comparing solutions of the 
stochastic equation to solutions of the corresponding deterministic equation. □ 



This paper is primarily concerned with the study of a diffusive version of the 
LSW model and the convergence of its solutions on any finite time interval as the 
coefficient of diffusion e goes to to solutions of the classical LSW model. This is 
a first step towards showing that solutions of the diffusive LSW equation behave 
similarly to solutions of the LSW equation for arbitrarily large time. In particular, 
one expects that solutions of the diffusive LSW equation coarsen in a fashion similar 
to solutions of LSW. Thus the average cluster volume is expected to increase linearly 
in time at large times, and this should be uniform in e for small e. In Theorem 3.2 
we obtained a uniform time averaged upper bound on the rate of coarsening. 

Although the effect of the diffusion decreases with time, discussion in the physics 
literature suggests that it plays an important role in the asymptotic behavior of the 
solution by acting as a selection principle. Thus at moderate times the diffusion 
produces a Gaussian tail to initial data (with compact support for example). It is 
the tail which determines asymptotic behavior, singling out the infinitely differen- 
tiable self-similar solution of the LSW model. 

The paper begins with a derivation of the diffusive LSW model from the Becker- 
Doring model. The least justified step in the derivation is the imposition of a 
Dirichlet boundary condition, which cuts off interaction between clusters of large 
volume and clusters of 0(1) volume that are almost in equilibrium at large time. 
The boundary condition is imposed at a particular cluster volume. One can see that 
this is a significant simplification since it implies that after the condition is imposed, 
the monomer density is for all subsequent times strictly larger than its equilibrium 
value. Despite this issue, the derivation of the diffusive LSW model from the 
BD model does offer a strategy for attempting to understand the mechanism of 
coarsening in the BD model. In order to carry it through one will have to understand 
in a precise way how a boundary condition is imposed on the diffusive LSW model 
by the BD dynamics. 

There has been some previous literature showing in an almost mathematically 
rigorous way a connection between solutions of the LSW and BD models. The key 
assumption required in this work is a type of upper bound on the coarsening rate 
for the BD model. It is not clear how this is related to the boundary condition 
assumption discussed in the previous paragraph. One should note however that a 
positive lower bound on the monomer density minus its equilibrium value, which 
is a consequence of the Dirichlet boundary condition, is a type of upper bound on 
the coarsening rate. 




6. Discussion 
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